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SUMMARY 

In the present report the theory of free turbulence 
propagation and the "boundary layer 'theory are developed 
for a plane— par allel free stream of a compressible fluid. 
In constructing the theory use was made of the turbulence 
hypothesis by Taylor (transport of vorticity) which gives 
best agreement v^ith test results for -problems involving 
heat transfer in free jets. 

The theory developed here considers two kinds of floi/: 

1, The boundary layer of a jet with temperature 
different from that of the surroundings and velocities 
that are small by comparison with the velocity of- sound 
(Bairstow number Ba <: 0,5), 

2, The boundary layer of a jet of high velocity and 
at a temperature eq.ual to that of the surroundings. 

The first deals with compressibility effect arising 
from the difference in temperatures inside an"d outside the 
jet; the second with the compressibility effect arising from 
the high flow velocities. 

The results indicated that the compressibility had' 
only a slight effect on the properties of the free jet. 
Furthermore it was found that a drop in the. jet., temperature 
had approximately the same effect on the properties of the 
jet regardless of whether the reducbion was due to ar-ti— 
ficial cooling of the jet or to the conversion of thermal 
energy to kinetic at -high flow velocities. - 

Compressibility factors (r and S) are introduced with 
the aid of v/hich it is possible to reduce the -v-ariations in 
all fundamental properties of the" boundary layer under the 
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influence of compressi'bility ■ to simple" linear relations. 
The results ohtained are valid' "for flow velocities up to 
the velocity of sound and for considerable temperature 
differences (up to +100 - 150° C) , 

lUTEODUCTlON 

In 1926 Toilmien*s paper was pulDlished (reference l) 
in X'/hich the author, on the "basis of the s eini— empirical 
general turbulence theory of Prandtl, developed the theory 
of so— called free turbulence — that is, turbulence in a 
free, hea,ted stream. In the same paper, making use of his 
proposed theory, Tollmien solved three problems on the 
propagation of free, heated jets: 

(a) - The boundary layer of an infinite plane— parallel 

Jet; 

(b) Plane— parallel jet escaping from a very narrow 

opening; 

■(c) Axially symmetric jet escaping from very narrow 
opening. 

About 3 or 4 years later (1929—1930) Swain's paper 
(reference .2) and Schl icht ing ' s paper (reference 3) extend- 
ing the theory of free turbulence to the case of the wake 
behind a body and developing the laws of flow in axially 
symmetric and plane wakes appeared. These laws are appli- 
cable to flows not too near the body, 

!Phe work of these authors was supplemented by experi- 
mental investigations of the velocity fields of flow. The 
use of one empirical constant enabled the above— ment ioned 
theories to be brought into excellent agreement with test 
results. In fact this agreement determined the success 
of the Prandtl— Tollmien free— turbulence theory and assured 
it v/ide theoretical and practical application. In 1935 
the article by Kuethe (reference 4) appeared in which an 
approximate method is worked out for the computation of the 
velocity profile in the initial part of a round jet. 

In 1935, 1936, and 1938 four papers by the present 
author were published in which Tollmien »s theory was ex- 
tended to the case of plane— parallel flow and axially 
symmetric jets escaping from openings of finite diameter 
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(an approximate theoryiof the initial part of the Jet was 
proposed); formulas were -worked out for the aerodynamic 
computation of the plane-^-pararlel jet; axially symmetric 
jet flow in the open .working part of a wind -tunnel with 
round and elliptic sections, hot and cold air Jets, and 
moreover, , methods v/ere proposed for comput ing> the air 
resistance of railway cars (in tunnels or on the ope& 
track), pipe systems and heat int e rchanger s (references 

6, 7), Ihese proposals and the flow theory itself 
were satisfactorily confirmed by test results and each 
year find wide application to engineering practice. 

The rapid development of the mechanics of turbulent 

flow has prompted the application of the physical model 
of the phenomena as conceived "by Prandtl and Tollmien to 
the solution of heat problems, those of the temperature 
distribution along the jet axis and over its c ros s— se c t i ons 
of heat diffusion from the jet to the surrounding space, 
and so forth. It is interesting to note that as a direct 
consequence of the Prandtl theory, in the case of a free 
jet and the wake behind a body, complete similarity is 
obtained between the temperature and velocity fields. In 
order to check this extremely important result Fage and 
Falkner (reference 8), in 1932, made measurements of the 
velocity and temperature fields in the wake behind a long 
cylinder of elliptic cross section. They showed that the 
theoretical velocity fields of Prandtl— Schlichting were 
very well confirmed by test results while there was no 
similarity of the velocity and temperature fields, and 
the heat transfer from the wake to the undisturbed flow 
is of greater intensity than follows from the Prandtl 
theory, ^ 

Taylor (reference 9) was the first to note the con- 
tradiction revealed in the free turbulence theory of Prandtl 
and presented a hypothesis according to which the tangential 
turbulent stresses in the flow were to be determined by the 
transverse transport of vorticity and :not by the momentum 
as proposed by Prandtl. The imperfection in the Prandtl 
theory was also pointed out in that it took no account of 
the local pressure gradients which have an appreciable 
effect on the momentum interchange but not on the, vorticity 
,tr,ansport. The ab.oye hypothesis, wit.h regard to the tur- 
bulence, was first proposed as far back as 1915 (reference 
10), On the basis of this hypothesis using only one em- 
pirical constant, as Prandtl did, Taylor obtained the vel- 
ocity and temperature profiles in the wake behind a long 
cylinder as exper imentaliy det ermined by Fage and Talkner^ 
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The Taylor theory of free turTjulence gave velocity pro- 
files accurately, the same as those given' by the Prandtl 
theory, and at the saae time removed the imperfection of 
the latter theory as regards application to heat problems. 
This. permits all the solutions of the problems in the field 
of flow meehailics that were based on the Prandtl theory to 
retain their validity. It madfe It necessary, however, to 
give preference to the Taylor theory for the further devel- 
opment of the problems of free flow and the wake behind a 
body, particularly in 'those cases when the problems are 
concerned with temperature profiles and heat transfer. 

In the present Report devoted to the further develop- 
ment of the theory of the free jet, a theory of free turbu- 
lence in a comprpssible gas is worked out and solutions are 
given of boundary layer problems of a free flow for the 
following two cases; 

(a) The boundary layer of a plane— parallel jet at 
small flow velocities with a temperature different from 
that of the surroundings, that iSj a nonisothermal layer; 

( b) Isothermal boundary layer at large (up to BaQ=l) 
flow velocities. 

It may be noted in conclusion that the free turbulence 
problems, in addition to being of interest in themselves, 
also possess a general interest since free turbulence repre- 
sents .the simplest case of turbulence free from the effect 
of viscosity. The study of free turbulence is a necessary 
preliminary stage in the study of turbulent flows in ^^eneral. 
It is, therefore, hoped by the author that the solution pro- 
posed in his present paper of the problems of free turbulence 
in a compressible gas possesses a certain usefulness for the 
study of turbulence" in other cases of compressible flow, 

I, EQUATION 01" MOTIOU FOR FREE TUEBULBirCB 



The problem will be restricted to two— dimens ional flow. 
In this case the differential equation of motion in the di- 
rection of the axis of abscissas assumes the f-ollowing form: 



bu 



bu 

+ -p-u — + pV 



5u 



P 



5t 



.by 



= - u- 



by" bx 



(1) 
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where 

u, V Inst ant aneous veiocity components 

p, P, \], instantaneous values of the density, pressure, 
and viscosity 

For the flow of a liquid of small viscosity alsout 
solid Tsodies the flow, as was shown "by Prandtl as far 
hack as 1904, may be divided into two regions; namely, 
a relatively thin layer of fluid lying close to the solid 
walls — the "boundary layer — in which the effect of the 
viscosity cannot "be neglected, however small its value 
may he, and the remaining part of the flow, in which the 
viscosity plays no part and' which . is therefore subject 
to the laws of flow of ideal fluids. The boundary layer, 
in turn, is assumed to consist of a very thin sublayer 
of purely laminar flow in direct contact with the v/all 
(no transverse turbulent fluctuations can be developed 
since they are dissipated by the wall) and a remaining 
turbulent portion of the boundary layer in which the 
effect of the viscosity may be neglected. Thus the 
study of the flow about solid bodies, the motion through 
pipes and, in general, of all fluid flows in the presence 
of rigid boundaries does not, in principle, permit neglect- 
ing entirely the effect of the viscosity. This circum-' 
stance constitutes the great obstacle in the development 
of the theory of turbulent flows. 

The distinguishing property of free turbulent Jets 
is the absence of rigid flow boundaries and hence of a 
laminar sublayer; This makes it possible to neglect 
entirely the effect of viscosity in all cases of free 
turbulence and explains the dynamic similarity of the 
jet flows — the npndependence on the Reynolds number — 
over a very wide range of Reynolds number. 

The differential equation of motion for two— dimens ional 
free turbulence may thus be written in the following form: 

5u , bu „ bu bP . . 

P bt * P^by =-bi 

Due to the quas i— s tat i onary state of the turbulent motion 
all its characteristics may be broken up into mean and 
fluctuating components: 



6- 



NAOA Technical" Kemorandnijn No-, -'1058 



•a. - u + u ' 
T = Y + yt 



P + P' 



P = P + P' 



(2a) 



On the average oreT certain finite tim-e interval, the 
fluctuating compLoneJit is. evidently equal.- to zero: 



•VI =■ u» ■= p » s P ' 0 



-(2b) 



In the -general caSe, however, this is hot -true of t"h'e \' 
squares 'of the f-luctuati ons and their ' products,' 

In eqiiation (2) th6 mean values and fluctuations are 
substituted for. th"S instantaneous -magnitudes. To average 
over the time, talce, account -of- conditions ( equati'on '(5b)) 
and neglect mom-ents of the third order-; - * . - 



« ■ • bu ' 



there iV obtained the differential eq.ua1;_ion for the average 
turbulent' flow of a compre s s ible fluid:' , .-. 



- Bu - - — bul f-- 7 , bu« - , bu' 

p u — + p V + p « u ' — + u p ' .+ o u ' 



ou' _ bu'l 



bP 

bx 



(3) 



To estimate the order of magnitude of the individual 
terms that enter the above equation it is not difficult to 
see that in the case of free turbulence all terms in th.e 
second bra-clc'ets ( 'vith' der iv-at ive s in"' "r e'sp'ec t tox) are -very 
small "by c oa-par is-on with the eorre'sp'oniing terms of the 
third brackets (iifith derivatives In -re-spect" to y) if •fche 
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■^r-incipal direction o.f ,the flow coincides with the axis 
of abscissas. Similarly one of the'terms of the third 
hraokets 

— buT 

V p' — 

^ bu 

is negligihly small. By neglecting the above small terms 
the following form of the differential equation of motion 
for the case of two— d imen s i onal purely turbulent flow is: 

pU-r-+p¥r— + pf V«-:-- + pV« - — = - ( 3a) 

bx by by by bx ' 



Free jets propagated 
liq.uid at rest, and wakes 
infinite undisturbed flow, 



in infinite space filled with 
behind a body surrounded by 
i, possess such small pressure 
gradients that they may be neglected. With this in mind 
the differential equation of motion for free turbulence 
in a compressible gas is: 



bu 

p u 



bu 

— + p T — + 
bx ^ by 



_ bu« 

p _ + p. 



bu"! 

byj 



(4) 



Ihe further steps in the solution of the problem depend 
on the choice of physical model for the turbulent flow. 
At the present time there are two models of interest in 
their application to free turbulence, namely, those of 
Prandtl and Taylor, With the aid of the Prandtl model 
results agreeing with experir.ent are obtained for problems 
in the field of jet mechanics (velocity fields, frictional 
stress, and so forth) but strong disagreement is obtained 
for heat problem solutions, (temperature field, heat trans- 
fer). The Taylor physical model gives the same solutions 
as the Prandtl model for the mechanical problems and further- 
more leads to solutions of the heat problems that are in 
good agreement with experiments. In what follows, there- 
fore, use is made of the Taylor model," ThV 'latfer is based 
on the assumption that the turbulent tangential stresses in 
the flow arise from the transverse transport of vorticity, 
that is, from the correlation between the vortex fluctuations 
and the transverse velocity components. In two— d imens i ona.1 
flow directed along the axis of abscissas the vorticity is 
given by 
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■ bu _ bV"] 
.by bxj 



(4a) 



where the magnitude 
negligihly small Toy 



of bT/bx in f ree_tur'bulence is 
comparison with bu/by so that 



1 bu 



2 by 



(4b) 



In its transverse transport, immediately before the loss 

of its individuality, the particle encounters a layer where 

the value of the vorticity differs from that in the layer 
from which it arrived by the amount 



by 



1 ^ h\ 

2 by= 



where is the mean free path of the fluid particle 

in the turbulent flow. The loss of individuality of the 
fltiid particle should be accompanied by a discontinuous 
change (fluctuation) of vorticity of amount 

o) ' = Aio 

Prom equation (4b) it is clear, however, that 

, 1 but 
2 by 

hence 



but 

"bT 



bfu 
bys 



(4c) 



With the. loss of individuality of a given particle there 
are naturally associated fluctuations of the flow velocity: 



bu 
by 



(4d) 



.19'ACA technical Keinorandum No, 1058 



9 



and of 'the fluid density; 

pt = It 



by . 



(4e) 



By making use of equations (4c) and (4e) equation (4) was 
reduced to the form 



or 



au — - bu 

pu — + p V — 
by 



— ou 
pu — + 



+ v It 1 



bu ^ ■ 

p T — + Y I 

by 



L. by2 



by byJ 



_<b_ r_ bu" 

by L Sy. 



= 0 



(5) 



(6) 



which is the differential equation o'btained on the basis 
•of the Taylor turbulence model. 



With respect to the nagnitude V' Ij it is necessary 

to make some assumptions hy which it is -associated with the 
velocity of motion and with the coordinates of the system. 
It is possible, for example, to make use of the generally 
accepted idea of Prandtl, namely, that the transverse 
velocity fluctuations are of the same order of magnitude 
as the longitudinal fluctuations: 



that 'is, 



7' ~ u« 



VI ~ Im 

- by 



(6a) 



Including the proportionality constant i.n the magnitude of 
the free path of the particle — the mixing, length (ij) — 
(equation (6)) is reduced to t,he more , simple, form* 



bu 
b^ 



bu 
by 



bu _b_ 
by by 



(7) 



*With. a view toward simplicity of notation the averaging bars 
over the letters are omitted in what follows so that p, u. 



and 



u. 



lip _are to stand for their mean values in time (p, 
and It). 
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In the special case of an incompres s i"ble fluid (p = 
constant) the following is ohtained 

bu „ bu . 3 bu b u . . 

bx by by by 

The corresponding eq.uation derived from the tur"bulence 
model of Prandtl for the free turhulence in an incompre&s— 
i'ble fluid was obtained by Tcllmien in 1926 in the follow- 
ing form (reference 1)} 

bu „ bu „ o bu b^u 
bx by by by^ 

To take into account the fact that the constant of pro-- 
portionality is determined from experii&ental data, it is 
seen that in the case of free turbulence the Prandtl and 
Taylor models give rise to the same eq.uation of motion. 

The value of the mixing lengths, as given by Prandtl 
and Taylor respectively, differ by the constant magnitude 



For the purpose of retaining the form of computation adopted 
by Tollmien and others the prandtl value of the mixing 
length is assumed. The differential equation of motion for 
free turbulence in a compressible gas then becomes 



Comparison of the above equation with the generally known 
equation for stationary flow 




(7c) 




(8) 



bu 

pu — - + pV 



bu 



&y 



bT 



by 



(8a) 



reveals the presence of "apparent" tangential stresses 
the magnitude of which is given by the equation 
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In the case of an incompressiblo fluid eq.uation {Q'b) re- 
duces to the, generally famili£5,r Prandtl law of turhulent 
friction: 



xy 



To solve equation ( 8) it is necessary to know a relation 

"between the mixing length \ and the coordinates of the 
system. The absence in the case of free flow, of rigid 
"boundaries that damp the fluctuating motions of the par- 
ticles led Prandtl to the assumption of constancy of the 
mixing length in the transverse direction of flow: 



l(y)= constant (9) 



It thus remains to estrihlish the law of variation of the 
mixing length along the axis of abscissas; 

I = l(x) 

The available experiment'al investigations of free flows make 
it possible without any particular difficulty to determine 
the form of the function l(x), A sufficient basis for this 
is' the experimentally established fact of similarity of the 
boundary layers in various cross sections of a given free 
flow (jets or wakes, reference 1), This similarity was 
revealed in a large numbor of experimental papers (Tollmien, 
Porthmann, Euden, Schlichting, and others) by constructing 
velocity profiles in nondimon s ional coordinates, for example, 
in the form of the relation 
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where 

u velocity g.t a point with ordinate y 

Uju velocity on the axis of the jet 

h width of the Jet (or of its boundary layer) -in tho 

given cross section 

The nondimens ional velocity profiles ( e quat i on (9a) ) were 
found to agree for the various cross sections. 

The similarity of the boundary layers at any two 
cross- sections of a given free flow must also he obtained 
with regard to geometric factors. In other words eq.uality 
is to be expected between the nondimensional mixing lengths 
for the various flow cross sections: 

-r~ - — «... = constant (9h) 



It is thus sufficient to establish the law of increase of 
width of Jet along the axis of abscissas in order that the 
law of increase of the mixing length be known. A very 
interesting consideration of Prandtl permits the solution 
of the relation b = b(x). It is shown by Prandtl (refer- 
ence 11) that the widening of the Jet (or of the boundary 
layer of the Jet) arises from> the transverse velocity 
fluctuations V' that is, 



dt 



( lOa) 



Because of the similarity of the velocity profiles at the 

different Jet cross— sect ions the following equation may be 
written 

5u ^ Um 

oy b 

and further 



db 
dt 



~ '^m 



(10b) 
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On the other hand^ the rate of expansion- of the jet 



dh 
dt 



dti dx 
dx dt 



(10c) 



that is, 



dh 
dt 



d;b 
dx 



Comparison of expressions (lOh) and (lOo) leads to the 
solution of the prohlem of the law of incre,ase in width 
of the free jet and of the mixing length in the flow 
direction: 



= constant 

dx 

h = X constant 



1 = c X 



(11) 



The Isjw obtained for the increase in the mixing length 
along the flow 



1 = c X 



(12a) 



is valid for free jets of various shapes; for the boundary 
layer of an infinite two— dimens ional flow, for a plane- 
parallel stream, axially symmetric stream and, in general, 
for those cases of free streams for which the flow profiles 
are similar. In the same manner, as previously described, 
the law of variation of the mixing length in a plane- 
parallel wake, axially symmetric wake, and so forth, may 
be obtained. In the present paper, which is devoted to 
the free jet only, no consideration will be given to wakes. 
By making use of the relation obtained for the mixing 
length along the jet the differential equation for free 
turbulence in a compressible fluid is reduced to the new 
form 



(12) 
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This is the general equation satisfying any case of a 



free jet 
the only 
"bulence, 
equation 
paper : 



of a compressi'ble gas. The magnitude c is 
empirical constant in the theory of free tur— 

In the special case of an inaompressihle fluid 
(12) assumes the familiar fpPDi given in Tollmien' 



&u 

u — - + V 
dx 



bu 
"by 



= 2c®x^ 



bu 



by by' 



(12b) 



II. DirFEEEKTIAL EQUATION OJ THE TURBULENT BOUITDAEY 
LAYER IN A PREB JET OP A COMPRESSIBLE GAS 
(PLANE- PARALLEL PROBLEM) 



Assume a plane— paral le 1 flow of compressi'ble gas 
extending to infinity in OY direction (fig. 1) with 



undisturbed velocity Uq , density Pq, and 

from the point 0, mix with the surrounding 
In accord with the law previously derived of 
increase in the width of the "boundary layer 
with the condition of similarity of velocity 
the velocity along any line Ocp drawn from 



start ing 
gas at rest, 
the linear 
h together 
pr of i le s 
the origin 

of coordinates 0 ( the latter coincided with the point 
where the boundary layer thickness h = 0) remains con- 
stant as will he shown. Prom the similarity of the 
velocity profiles it follows that the velocity at corre- 
sponding points of the flow are equal, that is, for 



h. 



5^3 _ ya _ 



• • • = constant 



there is 



u. 



^O ^O ^0 



: • • ' = cons tant 



But from equation (ll) 



h = X constant 



hence for 
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y/x ~ conjtant.. (iSa) 
there is the condition 

~ = constaiit (iSh) 



which' proTes what was required, since equation (l3a) is 
the equation of a straight line through 0, Thus in the 
turbulent boundary layer of a free flow the rays from 
a correspondingly chosen origin of coordinates are 
"isotachs, " 

The result obtained indicates that if the problem 
of the free plane boundary layer is solved in coordinates 
X and T\ = y/x the velocity will depend only on t\t 

u = UpfCn) (14) 



In order to eliminate the empirical constant from equation 
( 12) the following equation is set 

2cS = a'' (15a) 



and the following system of coordinates is chosen 

y 

X ; CD = — (15b) 
ax 

Then 

u = u^fCq) (15q) 

There is introduced, as is usually done for compressible 
flow cases, the stream function for the product of the 
mean velocity by the mean density and there is obtained . 

an ^ and p7 a — (l5d) 

cy bx 
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The density of the fluid in the case of free turljulence, 
for v/hich the pressure gradients may he neglected, de- 
pends only on the temperatures. 

The temperature fields in free flows, as shown hy 
the tests of Page and Palkner, Euden, Olsen, and others, 
are similar as the velocity fields. Otherwise expressed, 
the investigations of the temperatures show that in the 
free boundary layer the isotherms, like the isotachs are 
straight lines from the origin. Thus the temperatures, 
and hence also the densities, depend only on the non- 
dimensional coordinates ((p): 



e(cp) ) 

( ^ ( ^^^^ 
P = Po 1^^^) J 



On the hasis of the foregoing a certain function F(cp), 
the first derivative of which is equal to the principal 
component of the momentum, is introduced 



P ^ = Po ^0 ^' (17a) 

where and u^ are the density and velocity, respec- 

tively, in undisturhed flow. 

Hence 

\|/ = / pudy ^ Po Uq ax / r'dcp 
that is, the formula for the stream function is 

\i/ = axp^jU^F (I7h) 
and for the transverse momentum: 
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or 



(17c) 



Sulistit-ating expressions (l6), (l7a), and (l7c) in 
the differential equation (12) and making certain elementary 
transformations of the latter, using the expressions 



a = 3 



7 



2c2 : 



^ ax 



ocp 



x' 



btp _ 1 
l>y ~ ax 



(I7d) 



the differential equation of the turbulent "boundary layer 
in a plane— parallel free flow of a compressible gas in 

the form is 



ccp 



O CO 



o 
ccp 



K 



ocp 



(18) 



At all points of the bom^dsry layer except its cuter (cp^) 
ani inner (cp^) boundaries the derivative of the velocity 
is not equal to zero, that is, 

0 

c)cp 

Therefore the following is obtained 



ficp 



(19a) 



and further 



(iPb) 



Thus the differential equation of t.he turbulent boundary 
layer assumes the following form 
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To solve the above equation it is necessary to knov; the 
densiijy function k = K(cp), 

In the case of an incompressi'ble fluid for which 
K(cp) = constant = 1 and K« = ^K/bcp = 0 equation (20) 
reduces to the well— known Tollmien equation: 

F'" = I- (20a) 

corresponding to the case of the turbulent boundary layer 
in a free plane— parallel flow of an incompressible gas. 

In what follows equation (20) corresponding to the 
compressible gas will be solved for two cases: 

1, Tree flow at small velocities (up to Ba = 0,5) 

with a temperature differing from that of the surrounding 
space, that is, a noni s o the rmal flow. 

2. Isothermal flow at large velocities (up to Ba = l). 



Ill, BOUNDARY LAYEE OF IT OU ISOTHERMAL PLAFE-PAEALLBL JET' 
OP COMPRESSIBLE GAS AT l.,ODERATE ELOW VELOCITIES 
1, Heat Balance in the Turbulent Stream 



In the present section the flow of a compressible 
fluid (gas) in the boundary layer of a plane— paral 1 e 1 
stream at moderate velocities but at temperatures differ- 
ing froEi those of the surrounding fluid at rest shall be 
considered. In the later sections it will be shown that 
the effect of the compressibility of the gas arising from 
the high flow velocities is not large. Up to values of 
the Bairstow number of the order Ba - 0,5 — 0.6 the 
effect of the compressibility is barely appreciable, For 
this reason the equations and results which will be ob- 
tained in the present section devoted to the noni so thermal 
jet of small velocity will maintain their validity up to 
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velocities of the order of 0.5^ — 0.5 of the velocity 
of SOtLHd, T . - - 

To ohtain the law of temperature di s tr iljut ion in 
the hotindary layer of a free jet use is made of the 
differential equation of heat "balance, ;vifherein the 
molecular heat conduction .and the conversion of the 
energy of the viecous forces into heat is neglected 
with respect to the tur"bulent heat transfer in the 
same manner as the friction' due to the viscosity In the 
dynamic equation (2) was disregarded w.i.th respact to 
the turbulent friction. Then 



ST 5T &T 

p + pu + pV = 0 (21) 
ot dx ey 



where 



T temperature of the fluid 

t time, seconds 

It is convenient to "break up all the characteristics 
of the turbulent flow lui'.to their mean values sind fluctu- 
ations about the mean values: 

u=u+u», V = V+ V«, p = p+ p«, T=T+T» { 22a) 



so that on the average for a finite time interval the 
fluctuating components are reduced to zero: 



Vi = u» = p« = T'=:0 

Averaging with respect to time, while taking account of 
eqtiations (22a) and (22b) and neglecting moments of the 
third order: 



equation (21) is transformed into the differential equa- 
tion of heat balance for the turbulent quasi— stat ionary 
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f bT - hil f It _ ~"~iiT ^'""^^l 

p u-r' h p V--— + piu' r—+p u' -r — + up ' . 

■f-= ^' bf — otT' ' '■ ^o!"! 

+ ptV« + p V« + V p' ^ = 0 (22c) 

L • -by by • by J 

Assuming as in the case of equation' (4cL.) and density 
equation (4e) variations that the temperature change is 
the discontinuity at the instant of loss of individuality 
of the fluid particle transported "by the flow over a dis- 
tance equal to the mean value of the mixing length (ij) 
resulted in 

T' = Im ^ (22d) 
by 

Because of the presence in equation (22c) of fluctu- 
ations of the temperature gradients its further transform- 
ation hecomes impossihle. In order to eliminate this 
difficulty the equation of continuity of the flow is here 
resorted to: 

Fie + + rilS-+ ^^P ?') + b(p'u) ^' b(p'u')1 

Lbt bt J L bx bx hx bx J 

rb(pV) b(p7«) b(p'T) b(p«V«)-l 

+ + + i- + 1- = 0 (22e) 

I by by by by • J 

which, after averaging, assumes the following form 

bp b(pu) b(p'u') b(pV) b(p'V?) 

^ + , + ^ + — L_i + — L i. = c (22f) 

bt bx Sx by . by 

Subtracting the averaged equation of continuity from the 
instantaneous: 

^Bl + ^^P^') + b(p «u) ^ b(pY«) ^ b(p '?) _ ^ 
bt bx bx by by 

and multiplying "by T' gives the following: 
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bt bx bx by ^7 



whence 



ht hx bx by by 

bt bx bx by by 

By averaging the latter expression and taking into account 

the quasi-stationary state f -isi--! = p • -~ = the 

\ bt bt / 

following is obtained: 



pu 



dx 



\ djpu'T' ) , d(up'r) I I dj? VT 
1 dx • dx J'^L dv ~ 



T') , d{Vr/T') 



dx 



whence 



-- dT 



P» 



dx 



+ P V 



T, dT 



dy 



dx ~ 



L ^ 'd.v 



dx 



= 0. 



(23) 



Neglecting, in analogy to what was done with differential 
eq.uatlon (s), the small terns entering the second brackets 
and the term b(yp«T«)/5y in the third "brackets the differ- 
ential equation of heat balance in the turbulent flow is 
obtained in the sufficiently simple form: 



(24) 
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Jot the purpose of further converting this equation the 
relations shall be taken into account 

P=/t-^-. V=lj^, 7-=/t-^ (2d) 

ftnd the iars dropped from the notation, that Is, p = p, 
1S»u, T«!T, 7=V shall be set. Then 

dT . ,.dT, 



dx dy ^ 



(26) 



OX 



dT 



P dy dy dy dy + 



IS 



du dT 
dv dy 



(26b) 



The right hand side of the above equation gives the trans- 
verse gradient of .the turbulent heat transfer 

dWr . a dp du dl d \ , ^ du dT^ 
d\f dv' dy"dV~^~d^[^^'^-dj^-d^\- <27a) 



Trom this the expression for the heat transfer in the 
turbulent flow of a compressible gas is obtained: 

which, In the particular case of an ineompressible gap, 
bp/dy « 0 assumes the following foxm that 



, o du dT 

Wr ~ It' ^ -3 5 — 



(27c) 



In the case of free turbulence in a compressible gas, accord- 
ing to Taylor model, it is assumed 



(27d) 



the differential equation or heat balance Is written thus: 



dl 



= 2 c« 



' dx 
dp du dT 



dT 



dy dy d v dy 



(28) 
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For the turliulence model of Frandtl 



X = c X 



and therefore 

dT , ,.dT , ^ \ dp du dT , d / du dT\[ . 

Comparison of equations (SS) and (23a) shows that the 
Prandtl model gives a heat transfer half as great as 
that given by the Taylor model. Moreover, the Prandtl 
model, as shown hy Taylor, leads to similarity between 
the temperature and velocity fields; a result which is 
not dbtained from the Taylor model. In view of the fact 
that the results of Page's and Ruden's tests confirm 
Taylor's free turbulence model and refute the Prandtl 
model, equation (28) is used as a basis for this discussion, 

For an incompressible fluid (p = constant) the egua— 
tion of heat balance reduces to the following form: 



or . dT 



du dT 



dx ' ' dv d V \ dy dy 

2, Temperature and Density Distribution Laws 



(29) 



According to the results obtained in section ZI of 
this report 



p/t = PoUoF' (?); p V = p„ . «„ . {<fF' — F); p = p„x (cp); 

3 



_ y . 



ax 



3 r 

« = l/2c^ ; 



— F = 



/ F' 



d<f 



(30) 



It is assumed that the excess temperature fields (differ- 
ence between the temperatures of the stream and those of 
the surrounding fluid) in' the various cross— sections of 
the boundary layer of the free Jet are similar 



" ' 0 ' 0 ' atn 



(31) 
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where 
T 



^ata 
Then 



temperature along a raj <p drawn from the origin 
of coordinates (start of boundary layer of jet) 

temperature In the region of undisturbed flow 

temperature of the fluid at rest in the space 
surrounding the Jet 



AT ^Ti:^ AT AT '^^ 



(32a) 



Making use of expressions (sO) and (33a) to transform 
equation (28) it is found that 



and further 



.Fe' = y.' ^5 — - 

O<0 



EL 

X 



+ ■ 



d 



1^) 



H 

F' 



y. / H" 



The ahoTe equation on comparison with equation (30) 
leads to the differential equation 



(32c) 



(32d) 



"y. e' 



:0. 



(33) 



In the case of an incompressilile fluid (K - constant; 
K* s 0): 



In f»' = const; H' = const; 9 = Cjtp -|- c.,. 



(34) 



With e(«pi) * 1 at the inner boundary of the flow (in 
the region of constant velocity) and 6 (qpg) 0 at 
the outer boundary, the constants of integration are: 



whence 



^1 ("Pi - ¥^8) = 1; 
1 



NACA Technical hemorandum Ho, 1058 



25 



Irom this the law of temperature distriljution in the 
h^cundary layer cf a^ -plan.e^parall.el. flow, of .. aji , Inconv- 
pressihle fluid is found to be: 

cp ■- cpg 

e — ■ 

^1 - 9b 



The obtained linear law of temperature distribution 
is satisfactorily confirmed by Buden's tests (fig. 2, 
reference 12), 

A certain amount of disagreement with the tests 
occurs only at the boundaries where the temperature 
profile departs from a straight line and passes smoothly 
over to the boundary values of the temperatures. The 
therEial boundary layer is found to be somewhat thicker 
than the dynamic boundary layer. This fact is explained 
by the following reasoning. The linearity of the temper- 
ature lav; is obtained on the assumption of purely turbu- 
lent heat transfer with the molecular heat conduction 
neglected. This assumption was based on the analogy 
with the dynamic problem v;here the neglecting of the 
molecular viscosity led to a velocity profile which was 
excellently confirmed by tests, A specific character- 
istic of the velocity profile was that near the bounda- 
ries of the layer the velocity gradients and also the 
frictional stresses were so small that allov^ance for the 
viscosity could have no appreciable effect on the deform- 
ation of the velocity profile. In contrast to this the 
temperature field was obtained v;ith large gradients near 
the boundaries. This indicates that the molecular heat 
conduction at the limits of the dynamic boundary layer is 
of appreciable magnitude so that the temperature field 
departs from the straight line law and the thermal bound- 
ary layer will be thicker than the dynamic. Subsequently 
it is attempted to perfect the temperature distribution 
lav/ in the free jet by taking the molecular heat conduc- 
tion into account. For the computation, however, of the 
density and velocity fields in the boundary layer of a 
free jet such refinement of the temperature law is not 
justified since the accuracy of the density field will 
not thereby be apprecis/oXy increased while 'the mathemat- 
ical labor will be con,-, ic.era'.iy completed. 

On, the basis of the foregoing it was preferred 
to investigate the laws of flow in a compressible fluid 
without allowance for the effect of the molecular heat 



(35) 
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condTiction and to restrict the problem to the solution 
of the previous differential eq,uation; 



f.|:=0 <36a) 

whence 

In (k 8 « ) ss constant 

and 

K e ' = Cj^ (36b) 



The further -solution of .equation (361)) is predicated 
upon a relation between, the density and temperature 
functions. With this in miiid Clapeyron's equation is 
used J 

P - g R p T , ~^ ' 

) (-360) 
I'o =. S E Po To J 

which in the case of a free jet with constant pressure 
along and at right angles to the flow ( P = Pq = constant) 
leads to inverse proportionality between the absolute 
temperature and the densities: 

P _ 'T'o _ Tata + ATo . , 

(36d) 



Po T^ta + AT 



whence 



K = — 



(36e) 



Nond imens iona,l parameters, charact eriz'ing the degree of 
heating (or cooling) of the jet are introduced: 



Tata 



(37a) 
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where 

ATo excess temperature in the region of constant 

velocity (u = Uq ) 

ITg^^Q^ ahsolute temperature in the gas at rest surround- 
ing the jet 

This affords, in final form, the relation between the 
density and temperature functions in the "boundary layer 
of a plane— par allel turliulent flow: 

K = (37I3) 

1 + Gt 

Substitution of the ahove expression in differential eq.ua/- 
tion (36h) gives 



9 ' _ °i 



1 + te 1 + t 



(37c) 



or 



d(i + te) ^ , 

L = dtp (37d) 

1 + t3 1 + t 



Eq.uation (37d) is easily integrated: 

In (1 + te) = D^cp + XZQa.) 

The constants of integration (D^, D^) are determined from 
the boundary conditions given above: 

G(c?i) = 1 

e(cps) = 0 

resulting in cp-cp^ 

1 + te = (1 + t)^^-^« (38^) 
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and thus in the temperature distribution formula; 

AT (l+t)'^'"''^^-! 
0 = = i_- (39) 

The above expression substituted in equation (-SVb) gives 
the law of density distribution in the boundary layer: 

K = = (1+t) (403,) 



or 



« = _P, = (1 + t)'P^-^^ (40b) 



P 



0 



3, Development of the Differential Equation 



Consider the problem of the flow of non i s othermal 
jet at velocities that are small compared to the velocity 
of sound (up to Ba = 0.5), Por this condition the density 
profiles inay be considered practically indexiendent of the 
velocity profiles. The density will dexiend only on the 
temperature, the character of the dependence having been 
established in the foregoing as 



K = = (1 ^ t)^^-^- (41) 

Po 

The derivative of the density at a given point will thea 
be 



, . 9i-9s mil + t) 
B^=-- = -(l+t) — (42) 
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The ratio of the derivative of the dens ity . t o . the value 
-of the latter will he constant for a given value of the 
.Jet • temperature ' - 

k' ln(l •+ t) " " 

__ = ^ (^2^ 

Returning to the general differential eq.uation (20) of 
the houndary layer in a compressible' gas stream: 



jtii = _ p + ±. rUl. yil (44) 

its special form fo.r a noni s o thermal jet of moderate 
velocities is obtained as: 

p.n = _ p ^" i^iiJlil (45) 

After introduction of a special notation for the parameter 
which depends only on the temperature of the stream: 

S = (46a) 

The differential equation will then have the following form 

. _ i-m = _ r _ si>(f (46b) 

The above equ.ation is a common linear differential equation 
of the third brder whose general integral is of the form 

P = C^e"^^^ C,e^-^ 4- C3e^-^ (46c) 



whe* 6 



Gi » Cg, Cg constants of integration 

, kg, kg roots of the characteristic equation 
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+ ak*^ + 1 = 0 (46d) 



which in this instance reduces to the equation hy Cardan, 
As is usual for jet Isoundary layers equation (46c) has 
five lioundary conditions. 

1, At the inner "boundary of the layer where cp = 

(a) The gradient of the momentum is equal to 
zero ; 

=0 - that is, F"(cpi) = 0 ■ (46^) 

hp 

(h) The momentum is equal to the momentum of 
the undisturbed flows 

pu = pqUq— that is, J'Ccpj^) = 1 ^^Sg) 

(c) The transverse component of the velocity 

vanishes : 

Po^^o = 0 - ^^^^ ^(^i) = ^1 (463) 

2, On the outer boundary of the layer where © = cpg 

(d) The velocity gradient is equal to zero: 

^£^2 = 0 - that is, ^"(cps) = 0 (464) 
Sep 

(e) The velocity is equal to zero: 

pU = 0 - that is, P'Ccps) = 0 ^46^) 

The five conditions {4:6j__=,), are used to ascertain the thi-ee 
constants of integration 0^^ , O3, O3 and the values of the 
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nond imens ional coordinates of the-outer and inner limits 

of the t'oundary layer cp^ and cpg, !Po each value of the 

compressibility parameter S there corresponds certain 

values of the constants of Integration and the nondimen- 
slonal coordinates. 



4, Integration of the Differential Eq.uation of 

the Nonlsothermal Jet 

According to the foregoing the boundary layer of the 
nonlsothermal jet is characterized by the differential 
equation 

P" I + SP" + P = 0 (47) 
the integral of which is 

F = C.e^i^ + C^e^^^ + G^e^^'^ (48) 

The values of kj^ , kg, kg entering this integral are the 
roots of the characteristic equation 

k^ + Sk?' + 1=0 (49a) 
By means of the substitution 

k = — — (49b) 
the given cubic equation is reduced to the Cardan solution 



in which 



X* + 3PX + 2q = 0 (49c) 



P = - 

q = + ^ (49d) 



The roots of this equation are determined by the Cardan 
formula: 
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= n + V, X.a ^'-v^u + w^v, ' X^j --. jg^u + w^^v • (49e) 



the terms u land v ' "being given lay 



The coefficients ' and Wg are 'th'e conjugate imaginary 

cuTao roots of unity; 



1 . J's 



1 _ i 

2 2 



(49f ) 



la the present case: 



u = 



L 2 j y L. 2 J 



- / 



2 J 



- S' 



(oOa) 



Since u and v are real numbers, is the real root of 

the Cardan equation (49c) and and X3 are the conjugate 

imaginary roots. Correspondingly k^. is the real root of 
the characteristic equation C49a) and and kg the con- 

jugate imaginary roots. 



Setting: 

^1 = ai. kg = + pgi, k3 = - p^i 
inte.oiral (48) is transformed into 

H(p) = Cie^-i^ + Cge^°'5+Psi)cp + c 36^ ""fc-P 



(50b) 



(50c) 



As is known a pair of conjugate imaginary solutions of a 

linear differential equation of the third order defines a 

pair of real solutions expressed in terms of trigonometric 
f un c t i 0 n s : 



\ ^ ■ ^"^r y (50d) 
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In its final form the differential eq.uation of the bounds 
ary layer may be given in the following form: 



F (?) = C,e'' f + C,e'= ' cos + Cae"''^ sin (?,cp) 



(51a) 



The magnitudes a^^ , a,,, and are readily expressed in 

terms of the "thermal compressibility factor" of the Jetl 

5_ In (1+0 
9i — 'fs 

?or thl^ purpose the expressions (49b) and (50a, b) are 
resorted to, while taking into account the fact that S 
is a small magnitude of the order of 0.05 — 0.15. Than 



„=^_p., + f] + [s-f]>0; 



w, = 



27 

— 1— n'3 



whence 



A'l = « -f- — 

= WM -f- w^v = 



3+ 



2S-> 



27 

l + i|/3 
2 

1 — i VT 



2S'' 



27 
25-^ 



and further 



— ^i—^ 



1 — 



81 ' 



-~" 3 -[Y T+ 81 1 + ^^ 

3 ' 81 I * 2 



••' ~ 3 - Y 



1 + 
1 + 



25=^ 
81 

2S8 



81 



(51b) 



(51c) 
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A comjparison of equations (51a) and (51b) yields the 
formulas for computing , a^, and fig from equation 
(50e) for giren values of the compressibility factori 




(5ld) 



In the particular case of an Incompressi'ble fluid. 



- 0, it is 



«! = — 1 , «2 = 



2 ' 



2 



which is in complete agreement with the corresponding 
Tollmien solution: 



F==e- • + C.*- ^ cos [-'^'■pj +■ sin 



ri/ 



V3 



(51e) 



Zstlmating the order of magnitude of the individual terms 
in equation (51d) it is readily apparent that the terms 
containing ( S') can be neglected since S is usually 
considerably less than unity. Thus, for example, if the 
temperature of the jet Is 100°C higher than the surrounding 
temperature, S will have the value 0.1 and = 0.001. 

Thus without appreciable impairment of the accuracy 



3 



S_ 
3 ' 



3 -11 

1*2— 2 



(52a) 



The basic function of the boundary layer then assumes 
the following form: 



/=■(?) = e ^ ' I^C.e-" -\-Cy cos 



V'3 



+ C^e^ sin 



(52b) 



The expression in braces corresponds exactly to the Tollmien 
solution for an incompressible fluid: 



/'o(?) = Cie-M-Q« cos 



K3" 



4" sin 



(52c) 



The first derivative of T(cp), equal to the noBdimensional 
momentum, is given by 
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pa 



3 ^ 



— -^(Csl/3 — Cs)^^ sin 



[ 



,e sin 



^11 



(52d) 



The second derlvatire of '(cp), which is the nondimisnsioual 
velocity gradient, is 

F" (q,) = e" [ ' + 4" (Cs V'S - ^"^ cos [-"^ <p - 

+ Q 1/3) ^^^cos [J^ <p] - ^ (C, 1/3 - C3) e"^ sin [-^- ? J + 

" [c,^ c/^'cos ?] +C3e"^ Sin [J^ ^] ] (52e) 

For these considerations the factor ,S*/9 in the fore- 
going equation are neglected. 

There remains, on the basis of the five boundary 
conditions (46i_5), the determination of the constants 

of integration Cj^ , Cg, Cg and the values of the non- 
dimensional coordinates of the outer and inner limits 
(«Pi) ^^Pa) °^ boundary layer: 

= 0. -0 

The proTilem of predicting the basic function 7(cp) is 
then solTed. 

The five boundary conditions lead to fire transcen- 
dental equations with five unknowns (^i* ^21 ^3 » *Pi » 
cpa) solvable for any particular value of the compressi- 
bility factor. By applying the transformation of variables 
proposed by Tollmlen 

9 = 9 — 91, (53a) 

these equations can be considerably sinplified. The sub- 
stitution Is made so that 
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(53b) 



Then 



F(:f)=-e 3 I D,e • + D, e ^ cos cpj -f D, e ^ 
«nd hence, while hearing in nind that. 



Kin 

2 * 



(53c) 



9i = 0,,9, = ?,-?„/'o = ^i^' -r^D-.f' cos 
five equations with flye unknowns: 



y 3 



\ , "2 IV3-' 



^ 01 g '01 — 



01 



__ n. 
3 01 — 



■^oa g" ^ 01 — 0; 



which are slnplified to 



F — 



p " F' n 

' 1)2 Q ' 02 '^i 



2S 



(54a) 



01 - ri. /=■«,' ^1 + -^-; F^^"^-^; /^-o/r^ — 0,135; y='„," 0. (54b) 



With the first three equations of the above system the 
coefficients , Og, and of equation (53c) can be 

expressed in terms of qpj^: 



■A + ^(A+AV'3)=i+-|- 



A+4-(Ar3-a)=-^. 



(54c) 
(54d) 



from which 

9i-l 



' 1.5 9' "1^3+1.7 ^^^^ 



With equations (54c) the five equations of (54h) (with 
account taken of equations (53) and (52e)) can be reduced 
to a system of two equations with two unknowns: 



-D,e '■' + -2-(D, + A1^3)e 2 cos[i^(9,-cp,) 

-2[D^V2,-D,]e 2 sin|^i^(cp,-cpi)J=-0,13 5; (54f> 



ro, 



]- 
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~-^[Dys + D,] e^sin ^P, .(54g) 

where Dj^, Sj^, -and-Os are taken from equations (54c). 
Addition and su'btraction of equations (54f) and (54g) 
gire two cew equations of a somewhat simpler form: 

AK3 e^~cos [1^?(?,_=po] - 



_• ¥3 - ¥i ri/^'* 1 



2D,e — ^ cos 



- e 2 sin = 0, 1 3 5 



(55) 



Hext it is atteapted to determine the functional relation 
1)etween the constant* , D^,, Dji , cp^ , and tpg and the 
coapressihllity factor S making use of the fact, as 
will he shown later, that the compressibility of the air 
is of only slight effect on the free Jet. Putting 

A = Ao + AD,; Z), = D,- + AA; D, = D,, + AD,; | 

where Sio, Dso* ^sot <P'to« And cp^^ are known vi 
of the constants for the Ineompressihle Jett 



(56a) 



^10 = -212^ = - 0.0062, Ao = 2io±0:! = o,987, 



== 0,578, (pio = 0,98 1 ; 'f = — 2,04, 



(56b) 



and the saall increments (compare equations (86a} and 
'86c) with (84b))} 
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9 ^1.7 

>(56c) 

9 ^ 

Reverting to eq.uatlons (55) expanding into series while 

neglecting all terms of higher degree than the second and 
the products of the small increments and using the values 



e =4,51 e =2 0,4 

sin(-2.62) = - 0.5, cos(-2.62)= - 0,865 

results in 

0. 082S.= 0.380 Acpg - 0.255Acpi 
4,634S = 0. 130Acpg - IS.GOAcpi 

The solution of these equations leads to a functional 
relation hetween the deformation of the houndary layer 
and the c ompr e s s ihility factor: 

Acp-j^ = - 0,34S; Acpg = 0 (57a) 
which yield the corrections for the integration constants: 

AD^ = 0; ADg = |; AD3 = (57h) 
^ * 9 1,7 

The constants of the auxiliary function 5'(«p) are then 
equal to 

= - 0,0062, Dg = 0.987 + O.llS; D3 = 0.587 + 0.59S (57c) 

and the ordinates of the outer and inner limits of the 
"boundary layer 

9i = 0.98 - 0.34S, 92 = - 2.04 (58) 
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Equat ions •( 576 ) and (58) give the integration constants 
of the fundamental function iF(cp), The first three 
houndlary conditions of the given problem are used: 



■^(«Pi) = «Pi. S-tCcp^) = 1, P"(cpJ = 0' (59a) 



After suhstitut-ion of thevalues"^ P « and F" from 
equations (53"b), (52d), and (63^) the first ■boundary- 
condition gives 



^.1 



= e 



"5^1 |0ie^i+ Oge"2^cos [-^i] + Cse^^sin j^jZlcp 59 b) 



The constants of the compressible gas are again expressed 
in the form 



Ci = C^o + 4Ci, 0^= Ogo + ACg, O3 = C30 + AC3 (59c) 



vfhor e 



= - 0; 0176, = 0. 1337, C30 = 0.6876 (59d) 



are the corresponding values of the constants obtained 
by Tollmien for the particular case of an incompressible 
flow. In the same manner, according to (58) the coordinate 
of tJie inner boundary of the layer may he written: 



«Pl = 9i 0 0.34S , - (59e) 

whore cp^o ~ 0,981 is the coordinate of the inner boundary 
for the incompressible fluid. Substitution of (59c) and 
(59e) in (59b) and use of the series up to- the second power 
term (due to the smallness of, S) of thef expone-ntial and 
trigonometric expressions gives 
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e = 1 _ S^^^; ^-^i^e-^io (; 1+ 0.345]; 



= e 2 



[1 - 0.17S] 



cos 



s m 



cos 



l-V^io + 0.3S sin |^-:/_-cp^oj 



> (59f) 



— 0, 3S cos 



[-^10] 



The products and squares of small terms are disregarded 
and the following relation is taken from the "boundary 
conditions for the incompre s s ihle gas; 



«Pi 0 



whence the equation connecting the increase in the con- 
stants due to the compressi'oility with the compressibility 
factor St 



0,375 ACj^ + 1. OSACg + 1. 235 AC3 = 0.340 S 



(60) 



The second boundary condition F'(cpj^) = .1 in combination 
with equations (52d) and (59b) gives 



1=8 



S 



^-Cie'^i+ I (C3+ C3 >/r)e * COB -"^cpi] 



CP: 



sm 



-^1 



(61) 



Substitution of expressions (59c), (59e), and (59) in 
equation (61a), while neglecting products of small 
quantities, and making use of the particular form of 
this equation obtained for incompressible gas yields; 
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a second equation for the relation between the increments 
in the constants of the function' P(cp) and the compress— 
ilDility factor: ' ' ' 

-O.375A0i - 0,52ACg + l.SSACg = •0.654 S (62) 



The third houndary condition F"(cpj^) = 0, together with 
eq[uations (52e) and (61), gives 



S ^1 

0 = e 



"^^|-C,e (03^5- C,)e^ cos [^|^,] 



1 
2 



(C3+ 02^)6 « sin |^:|cp,J [l^FiJ (63) 



By substitution of expressions (59g), (59e), and (59f) 
in equation (63), while neglecting the products and 
squares of small terms and taking account of the fact 
that in the case of incompressible gas equation (63) 
asstimes the form: 



9io 



0 = C^^oe 



'P-.l (C3,X- C,,)e'= cos [4.o] , 



the third relation between the 
of the function S'(q))- and the 

0.375 AO^ - 1.6 ACg + 



increments of the constants 
compressibility factor is: 

0.32 AC3 = 0.332 S (64) 
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The solution of the system of three simultaneous eq,ua— 
tions: 



0.375 ACj^ - 1.60 ACg + 0.32 AC^ = 0,372 S, 



-0.3.75 AC, - 0.52 ACg + 1.55 AC, 



0.654 S 



0.375 ACj^ + 1.08 ACg + 1,225 AO^ = 0.340 S 



(65) 



gives the laws of variation of the- c-onstants of -integration 
of equation (521)) under the effect of compressibility: 

tC^ = 0, ACg =^-0.14 S, AC3 = 0.385 5 (65a) 



hence the integration constants: 



Ct_ = - 0.0176; Og = 0.1337 ~ 0.140 S; 



O3 = 0.6776 +■ 0.385 S 



(66) 



Substitution of expressions (SSh) and (52e) in (66) while 
neglecting the products of small quantities, the function 
P and its derivatives are obtained in the final form 

FsPp+AF; = + AF» ; 2"' = pi'^ + AP" (67) 



where Pp, I'd ' « -^o" '^^^ -values of the functions 

and its derivatives for incompressible gas ' ( Tollmi en ' s 
s olut ion) : 

qvw, - w. . -.-.w.,- '^cp + 0.6 876e^sin ^g-cp 



Po'(cp)= 0.0176e "^+0. eeSSe^cos 



9 



— cp + 0.228 e2 sin 



1^1 



>(67 



^^"i^) = -0,0l7&e ^+b.528e"^CDs 



: _ cp 

- 0.93 06^ sin ^cdI 
u 2 ^ 2 ■ J 



and AP.-AP', and A?" are the. increments of the function 
and its derivatives under the influenpe of compressibility. 
The increment of the function is 
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- r /" 1 ^ Jz 



(67b) 



that of the first derivative 



AS"' = 



1 

+ P_ ^ 0, 8 Oe^ cos 



1 

- 0, 93e^sin 



■cp 



|(67c) 



and that of the second derivative 



CP- 



AF" = 



S 

3 



cpT "+2P « - 1.2e^ cos 



V3 



-cp 



9 

+ 0. 22e^ sin 



V3 
2 ^ 



(67d) 



In conclusion it should be noted that the author's 
carefully conducted numerical 'solution, based directly 
on equations (55), showed couiplet,e 8?veement" with the 
functional results, (67a) aAd ( 6 7a , t\pp r ox imat e 1/ obtained, 
notwithstanding the fact that for the numerical solution 
a very large value was chosen for 
factor (S = 0.135) corresponding 
jet has a temperature 150°C above 
ature , 



the compressibility 
to the case where the 
the surrounding temper- 



It vfill be recalled that the first derivative of the 
function F is the momentum of the flow in the direction 
of the X— axis: 



Po^o 



= pt 



(68a) 



To obtain the velocity u/ Uq it is necessary to apply to 
the law of density distribution in the boundary layer: 



CD. 



9 



Po 



(68b) 



then 



u 
' 11 



(68c) 
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The momentum velocity ratios in the direction of the 
y— axis are 

^ = a (cp?' - F); = (6 3d) 



Po^o ^0 K 



To facilitate the computation of the components of 
the velocity and other characteristic magnitudes of the 
nonisothermal layer, tahle I is appended with the computed 
values : 

I'd. Fo', Fo". ^ I AF, I AF'; | AF" 
The tahles 2, 3, and 4 contain the values; 

F ' coF» ~ F F'2 
F. Ft, F", K, ~; 2——; cpF' - F ~, 9 

from which the velocity profiles, velocity heads, den- 
sities, and temperatures may Tae computed for the follow- 
ing values of the compr e s s ih il i ty factor: 

S = - 0.074; S = 0. 06 05; S = 0.1115; 



the given values of S corresponding to the jet temper- 
atures of —60°, 60°, and 120° C, respectively, above 
the surrounding temperature. 



5, Fundamental Properties of ;Boundary Layer 

of Honisothermal Je'' 
(a) Geometry of the Jev 

The nonisothermal jet, that is, having a temperature 
other than the surrounding space has, as explained in the 
foregoing, the interesting property that its outer Taound— 
avj (u = 0) remains constant for variation v/ithin very 
wide limits of the temperature increment (At = Ji 150° 0): 



cpg = constant = — 2. 04 (69) 
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The inn.ep. ■boundary.- ( u <=/ vlq) expand-s ■ somewhat when the 
jet is cooled and compresses v;hen heated: 

cpi = 0,981 - 0.34 S (70) 



vfher e 

ln( 1 + t) 

S = compressihality factor of the jet 

t = Atp/Tata ratio of the temperature increment (At in ° 

to absolute t emperature ■ ( T^^a. ~ 273 + 
*ata° surrounding fluid. 

According to formula (70) the folloiiring results are 
obtained: At = -6 0°; cp^^ = 1. 005; At = + 6 0°; cp^ = 0.960; 
At = + 130°; cpi = 0.942; At = 0° - cp^^ = 0.981.- 

The nond imensi onal width of the jet boundary layer 
depends on the temperature difference above the surround- 
ing temperature in fche following manner: 

b= — = 3. 02 - 0. 34 S (71) 
ax 

The boundary that separates the initial mass of the 
jet from the entrained mass is determined, as is known, 
from the condition that at the partition surface (cpg) 
the stream function = 0, or, what amounfs to the same 

thing, I'(cpjj) = 0. The relationship between the boundary 
of the core of constant mass flow of the jet and the com- 
pressibility factor v/as obtained as follows: 

In the case of i'ncorapre s s i ble fluid it is 

^.(cp^^) = 0 

For a compressible fluid 

J' = + AF. ^. = cp^^ + Acp^ 

then 

i'(cp^) = ^^(cp^) + A5'(cp^) = Fo^930 + ^^3^ ^ ^^^^30 + ^^3) = ^0(^3^) 
+ Acp^Po'(cp^^) + Ar(cp^^ ) + Acp^AF«(cp^^) = 0 
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Omission of the small magnitudes of the second order 
(Acp, AF') while ."bearing in mind that 

F„(cp ) = 0 



0 '30 



gives 



hence 



S 

AF(cp J 3 
^jp = ^z.a.c.l_ s S 0. 29 S 

° 30 

The result is the formula for the boundary of the constant 
mass flow core: 

cp = - 0. 185 + 0.29 S (73) 

3 

The same result is obtained by direct interpolation from 
tables 2, 3, and 4, 

(b) Velocity, Temperature, and Density Profiles 

In the representation of the velocity profiles the 
magnitude 

- (P~^s 
cp = 



serves as ordinate, the effects of heating and cooling 
of the jet can be compared from the plotted velocity 
distribution curves. The curves shov;n in figure 3 are 
for At = + 60° C, -60° C, and + 120° C, 

On figures 4 and 5 are shov.'n the density and temper- 
ature fields in the Jet boundary layer for At = — 60 
and + 6 0° C. 
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Figure 6 gives the field of veloci ty ^heads (P'^/k) 
for At = - 60° C; 

p n. s 



fields for At = - 6 0° C: 



Figure 7 finally gi'Ves the fields of transverse velocity 

0° C: 

V cpF ' - F 



The corresponding fields ohtained "by TollTpi'^.-a for the 
iiac orjpr e s s i b 1 e fliud (.^^^q ~ '-') sho^'ii in ''iguies 3 

to 7 for comparicciia (Cee tahles II, III, and IV.) 

(c) Rate of tiass Flow 

The quantity discbF.vged per second in the turhulent 
c ouipr e s s ihl e jet d.3s not ei-C'rely correspond to the valu.es 
of the stream function on acooant of the fluctuations in 
the density. The mass flow per second is 

m = fp udy = fp udy +/ p'u'dy=\l/+ /p'u'dy; (73) 

But the stream function is determined by the expression 

^ = axp^u.pF 

while 

/F 

— — - _ hp bu Sk; \iC J 

SO that the mass flow pe-r second is 

'F' 



m = axpoUo F + a / « i 

Sep 



F + a / k' — ^i-^— dtp, + constantj (74) 
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or in nondiinensional form: 

, €) 

, I 



J = = p + a / n» — dcp + constant (75) 



In particular the nondimensional value of the entrained 
mass of fluid sucked into the free jet from the surround- 
ing space is* 




(76) 



The nondimensional magnitude of the retarded mass of 
fltiid in the initial jet core is 



3 

n 



+ a. J y:^ 



9cp 



(77) 



The values of 1^ and may be olJtained from tables of 

integrals and "by approximate integration (for example, by 
Simpson's method) the v&lues of the integrals: 



A, 



— dcp 



(73) 



for various values of the temperature of the fl0'^r, 5'urther- 
more, assuming a few values of the turbulence factor a, 
for example, taking a = 0.0845, according to the tests of 
Tollmien, it is not difficult to determine the values of 
the mass of the initial jet and of the entrained mass for 
various values of the temperature, 

*The subscript i hereinafter refers to the boundary of 
the region of constant velocities; subscript z to the 
region of fluid at rest; subscript 3 to the core of 
constant mass flow. 
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The integral in expressions (76) and (77) may "be 
integrated by parts: 



K 



dcp = — 
8cp « 



~ ?'cLcp 



3 3 

Prom the "boundary conditions it is known that 

P«a = 0 



Moreover, according to eq.uations (41) and (46a): 

ln(l+ t) 



K=(l+t)^--^^ S = 



CPl - CPs 



henc e 



(cpi-<p) S 
tt = e ; k: ' = 



(79) 



Omission, as in all parentheses of the preceding section, 
of the terms with factor leaves 



ocp 



dcp = - -- p« 
ft 



= - SP«, 



(BP) 



whence the expression for the entrained mass: 



^2=- C^-Ccpa) - aSF«(cp^)] 



(81) 



A glance at tahles II, III, and IV corresponding to 
particular cases of the noni s othermal Jet: 

S = _ 0. 074 (Atj3 = - 60° C) 
S = 0. 06 05 (Ato = + 6 0° C) 
S = 0.1115 (Ato = + 120° C) 
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shows the Talues of , corresponding to the 

partictilar values of the compressibility factor S 
on the assumption of a = 0.0845 for the coefficient 
of jet structure (turbulence) — according- to the data 
of Toilnlen and CAHI Therefore 



a S 7i(<p ) « 0. 05 S (82) 

The values of '(q>s) ?a * ~- 3. 04 is computed 

according to formulas (67) and (67b): 

r(q»a) ■ - (0.388 - 0.27 S) (83) 

Lastly equations (82) jsnd (83) added together give the 
final expression for the nondimenslonal value of the 
entrained mass of the Jets 



■a " 0.388 + 0.22 S (84) 
The Initial mass of the Jet is 



i — = /«.=/•. 4-a / y/ — — — 



x' |3 /• y" 

= r^a.^.F\-aj—.F.d, 
According to (79) ' 



(85) 



y X 
y. y. 

h'enee 

= + aS [F\ — F',] -Tr Fj — 0,034 5. (86a) 
Turther, equations (67) and (67c) give 

/=", =^ 0,981 + A/=-(^),)r>? 0,98 + 0.02 5. (86b) 

Substitution of (86b) in (86a) yields the final expression 
for the nondimenslonal mass of the initial part of the Jet: 

■ Wi-- 0,-981— 0,01,45. , (87) 
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The total (nondimensional) mass flow in the 'boundary layer 
is eaual to 



/«, + /«, = 1 .369 + 0,2 S. 



The smallness of the value of S Indicates that the 
mass flow per second for the non i sothprmal boundary 
layer differs very little from that for the Incom— 
presslhle Jet. 

(d) Jrictional Stress* 

In section I an expression was derived for the 
frictionf.l stresses in the turbulent boundary layer of 
a compressible gas, which according to equation ( 12a) 
can b e written 



= 2 c^x^- _f 1^ ^ (p i^yy + const. 



(89) 



But 



b'^nce 



ax 



— Y> 



d 



F' 



dp + const. 



The nondimensioual value oi tne frietional stress is 




F' 



£i?-|- const. 



(90) 



After corresponding transformations Dartial integration 
yields 



C)3 



J L ^9 



+ 



rf-f + const. 



(91) 



J he values of T of greatest interest i.s that at the 
boundary of the core of constant flow T = 7^ since it 
determined the energy loss due to suction of the entrained 
mass. This value of T is evidently given by 



F' 



(92) 
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At the inner limit of the boundary layer 
hence 



0, 



In addition 



where according to (79) 
whence 

F' 



— ) 



0^ (fi - 9) 

Since the values of S are small 



it results in 



—±-L = F" + ^(9 - P" + SF', 



J 



L d:i 

F 



= /=•"- -h 2 5 (<p — cpj) 4- 2 
= /•"•"••' -f 5" (9 - cp.) P"-' + 2 

1 



df = S J F^-d}. 
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Vlth the aid of the approxlaate lzit.egrat ion method and 
table I the Integral beeoaea 

4 

3 

within the limits ep > - 0.185 to qp,^ > .0.981 (tne 
girem Taluet of «_ and fp^ correipcnd to the incoa— 
preeeible Jet), koreoTer,at the point tf (for Tarioue 
ralvee. of S): 

whence 

~, = 0,27 — 5 0,32 -|- 0,62 54- 0,4 5. 

The final ezpreeeion for the nondlnensi onal value of the 
frletional streea at the boundary between the initial and 
entrained aaeeea of the Jet 

= --i— = 0,27 + 0,7 S. (94) 
"Po-2- 



(e) Heat Transfer 



In section 71 (par. 1} the differential equation 
was obtained for the transrerse heat transfer due to the 
turbulent fluctuations in the free Jet of the compressibl 
Sasi 



1 



dW. 



£?P du dT d 
dy ' dy ' dy ~^ dy 



dii dT\ 
~dy ~d^) 



(95) 



The heat transfer from the initial mass to the entrained 

1% studied 



■ass across 1 a' of partitidn surface ^ 
next: " 



«',=C,.^,2..=..[J^.^._ 



d'l , , du. 



dT 



3l 



(95a) 



Troa prerlous results 



d'i ~ 'o r~ 



du 



= K„-e<'f-''^^ [F"-\-SF'\, 
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whence 



3 



(c|>-cpi)S 



dcp 



(96) 



Use is also made of the known relation 

+1 



1 + t e 



After certain simplifications, while neglecting very- 
small terms, equation (96) assumes the forms 



Cp S Po ^0 AToa 



SP« 



9i 



+ S F" + I"' (cp- cpJS+SF' 

-93 J 



, 3 



Further transformations give 



(961)) 



aCp~s"po~^o^^o 



9i-92 9i-9s 



P'Kcp^) 



(97) 



But* 



2'"(93) = Fo"(93) + ^^"(93) = J'o"(93o) + ^9 32"" (930) 



+ AF" (cD ) = F„« (cp ) 

■ - O ^30 



30 



and, further, 



9i = 9io - 



9„ 



0. 185 + 0. 29S 



*In region cp the magnitudes F'" and AF'" assume 
values of the ''order of zero. 
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The foregoing relations yield 

— = Ji 120^ + o,34S- 

agOpPoUoATo cpiij. - (p^o 

where 

^o"^«P3o) = 0.52,- 0 = 0.981,,. cp3 o= -2.040 

Thus .the nondimensional magnitude characterizing the 
heat transfer through the boundary of ^ the Intial mass 
of the nonlso'thermal Jet is equal to - ' 

W = = 0.173 + 0.34S (99a) 

^ agCpp^UoAT^ 

The coefficient of' heat transfer from the initial mass 
to the entrained mass equal to the amount of heat trans- 
ferred per hour through 1 square meter per degree differ- 
ence in temperature is 

a„ =.-i p 26 00 a s C_, p. Uq [0.172+ 0.34S] — (99^3^ 

^ . ATq " hr °C 

For a turhulence coefficient a = 0.0845 and g = 9.31, 
it yields 

a,, = 3000 Cp pp UqC 0.178 + 0.34 S) (lOOa) 
For air with a specific heat of the order of 



it is 



Cp = 0.24 



C £1 1 

a.m = 720 po Uo (0.172 + 0.34 S) g- ( 100Tb) 

m2' hr C 



The effect of the compressibility on the heat transfer in 
the nonisothermal jet is appreciable. For " ejcample , for a 
temperature difference ATq = 60° C (S = 0,07) the heat 
transfer in the compressible jet differs from that .of 'the 
incompressible fluid by, about 15 percent. 
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(e) Concluding Hemarks on the Noni so thermal Jet 

As is seen by the previous discussion the effect 
of the compressibility of the fluid (gas) on the fundar- 
mental properties of the boundary layer of a nonisothermal 
flov/ is insignificant. In particular, on lowering the 
jet temperature 60° C below that of the surrounding 
medium (this corresponds to an increase in the velocity 
to Ba = 1. 0 in a jet of high velocities) the angle of 
divergence of the boundary layer increases by 0.7 per- 
cent; the angle of dissolution of the core of constant 
mass flow increases by 11 percent* the nondimensional 
value of the entrained mass by 3,7 percent; the mass 
of the initial jet remains practically constant; the 
frictional stresses decrease by 16,8 percent; the heat 
diffusion is reduced by 15 percent. 

The results obtwinea are evidence of the mainte- 
nance of the dynamic similarity of the jet for appre- 
ciable changes in its temperature,* It has further 
been shown that if S = ln(l + t)/(cp^ — cpg) is taken 
as the compressibility factor of the nonisothermal jet, & 
vrhere t = Atp/Ta,ta ratio of the temperature 

increment of the jet to the absolute temperature of 
the surrounding medium, the change in the fundamental 
properties of the jet with S is linear except for the 
nondimensional coordinate giving the dissolution of the 
outer boundary of the jet which remains unchanged. It 
will be shown later that the cooling of the jet at small 
velocities has the same effect on its fundamental prop- 
erties (velocity profile and friction) as an increase in 
the Bairstow number. 



IV. BOUNDARY LATER OF A PI.A1IE-PARALLEL JET 

AT LARGE VELOCITIES 
1, General Considerations 



The investigation so far involved the case of a 
free turbulent jet having a temperature different from 
"that of the -surroundings and a velocity small by com- 
parison with the velocity of sound; that is, the effect 

*The simil?.rity of the jet for a wide, var iat ion in the 
Reynolds Q-umber was discussed in a previous report (ref- 
erence 6), 
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of corupr essi'bility arising from the difference in. tem- 
perature alone was studied. The follo^ving deals with 
jets of large velocities (up to the velocity of sound) 
under the condition that the temperature in the reservoir 
from which the jet escapes is equal to that of the sur- 
roundings, or otherwise expressed, the effect of compress— 
ihility due to high flow velociti-feT will he investigated, 

' 2, Derivation of the Density Punction 

The air temperature ahead of the nozzle (in the 
region of small velocities) is equal to the temperature 

of the surrounding medium. In this case there will be 
no heat transfer between the jet and the surrounding 
space so that the heat content of the air will be uniquely 
associated with the flow velocities (the energies of the 
pulsating and transverse motions are neglected): 

Cp(r - TJ = ~-~ (uo^ - u^) (101) 



vrhere 

Cp specific heat at constant pressure 

A = -^T,-;^ heat equivalent of mechanical work 

g acceleration of gravity 

The relation 

give s 

_T_ 
T 

o 

The velocity of sound in the region of undisturbed flow 
read.s the value 



(lOl) expressed in nondimens ional form 



1 + °1_ 
2gCpTo 



u 

1 - ( — 
^u 




(101a) 



°o = /-H- g ^ ^0 



(lOlb) 
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where 

specific heat at constant volvime 
Cp/c^ = k adiabatic coefficient 
B gas constant 

The Bairstow number in the undisturbed region is 

Ba^, = Uo/co and AH = Cp - 

hence 



= 1 + 



— 3a/ll-(^_J j (102) 



The foregoing expression gives the relation between the 
teii;perature of the flow and the velocity. Particularly 
in the outer boundary of the flow (region of air at rest) 
where u = 0: 

^BIB: = 1 + Ba-* 

At the inner boundary of the boundary layer (undisturbed 
flov/ xtfith velocity Uq) : 

T 

A thermometer, however, mounted at any stationary point 
of the flow will show the same temperature — the 

temperatixre of the air at rest — since the velocity of 
the flow drops to zero directly at the wall of the ther- 
mometer. Thus a stationary thermometer in the flow shows 
not the actual temperature of ~TE"n] ow~"but the temperature 
of the retarded air, the stagnation temperature. 

The above flow considered from the point of view of 
the stagnation temperature is isothermal as a result of 
which it was possible to assume that no heat transfer 
exists and to apply the heat equation in form (lOl), 
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The foregoing prolDlems of free turbiilence involve 
isolDars (the pressure gradients are negligi'ble and not 
taken into acc oun t );. henc e the density is inversely pro- 
portional to the aissolute temperature: 




(102a) 



Therefore the density function is: 



K (cp) = 



k-1 s r 



2 



k(cp)/ 



(I02b) 



The solution of the above equation for K gives the 
density function in the form 



1+ Jl+ 2(k- l)BaQ^ 



/ k- 1' 
Ba^ ( 1 + ~g~Bao ' 2 



k- 1 ^ 3 



Ba 



0. 



(103) 



and the introduction of the special function of the 
Bairstow number: 



r = 2(k- 1) Ba, 



1 + 



whence the calculation of the derivative of the density 
function affords 

^1 rPip" 

/iTrpTi [ 1 + yrr^ ^^^^^ 

3, Derivation of the Fundamental Differential Equation 

With equation (104) the general differential equation 
(20) is reduced to the following special form which satisfies 
the problem in question. 
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= - '^ '^ 11- ] (105) 

Sep I J1+ rr»^;[l rF's] - . 

whence 

jtM-_j' + _£.ri''t 

&cp L yr + r F ' » 



and 



b r 1 

- P = (105a) 
&cp ! / 1 + rP ' s J ' 



Further transformations yield a differential equation 
of the form 



P m s= J- F yr+rF « ^ + ( 1 06 ) 

1 + rF'S 

In the special case of incompressible flow xirhen r = 0 
(equation (106)) reduces, as expected, to the known 
Tollmien equation: 

F'" = - F (106a) 



Differential equation ( 1 06 ) as well as the Tollmien equa- 
tion (l06a) contains the five "boundary conditions (46i_s) 

which are applied to obtain the three constants of inte- 
gration and the values of the nond imens i onal coordinates 
of the outer and inner limits cp^ and cpg of the boundary 
la^er. To each value of the compressibility parameter ( r) 
there corresponds the values of the integration constants 
and nondimen si onal coordinates. 

Since the functional solution of equation (106) is 
impossible, it is necessary to apply the method of 
numerical integration to each particular value of the 
compressibility parameter (r) or what amounts to the 
same thing, the Bairstow number ( Ba) , The most suitable 
method for solving the given equation appears to be the 
Adams method which gives good agreement for the given case. 
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4, Huj^erical Integration of Differential 
Equations "by the Adams Method 

If the values of a function and its derivatives 
at a certain point are known, the values at a neighbor- 
ing point (a + ,h) can he obtained with the aid of the 
Taylor series. After the values of the functions and 
its derivatives at the point (a + h) are determined, the 
values at other neighboring points (a + 2h) , and so forth, 
can then be found. In general ' if Yi^ = Y( a + kh) is known, 
then 



Thus, passing successively from point to point, it is 
possible to compute a table of values of the required 
integral, that is, of the required function Y(x) over 
the entire integration range. The smaller the size of 
the interval h the more accurately are the values Yjj 
determined although, on account of the large number of 
intervals, the accumulation of erroi's may become con- 
siderable unless a sufficient number of terms of the 
Taylor series is employed. The fundamental disadvantage 
of this method, which was proposed by Euler, is that it 
is necessary to compute the higher derivatives which, 
for arbitrary form of the function Y may yield very 

complicated expressions. In such cases a considerably 
less complicated method is that by Adams in which the 
increment in function at a certain interval is expressed 
by the first differences in the so— called linear incre- 
ments of the function, the second differences, in the 
neighboring intervals , and so forth. The Adams method 
does not require the higher derivatives and gives good 
accuracy even for large integration ranges. Particularly, 
as v/ill be shown later, in solving the differential equsr- 
tion of motion for the compressible gas jet (equation 
(l06)) — the Adams method gives very good agreement. 

The linear increment in a certain interval (a, a+h) 
is given by the product 

'In = h Y'^ (a) 
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The first difference of .th^ linear increments of the 
function in neighhoring intervals- is given "by 



n— 1 n .n— 1 
The second difference "by- 

'n— s 'n— I 'n— a 

the third difference "by 

'n— 3 n— g n— 3 

and so forth. 

In general if there, is a tatle of values of the 
required function at n + 1 points it is possible to 
draw up a table of linear increments at n + 1 points, 
first differences at n points, second differences at 
(n — l) points, third differences at (n — 2) points, 
and so forth. The Adams method makes it possible to 
compute the value of the function' at the (n + 2)th 
point, the value of the linear increment at the (n+ l)th 
point, the value of the first difference at the nth point, 
the value of the second difference at the (n — l)th point, 
and so forth. By the same method it is possible to proceed 
to the (n + 3)th point, and so on up to the end of the 
entire integration range. 

The extension of the integration table from one 
interval to the next by the Adams method i.s effected on 
the basis of the following considerations: 

1, The Taylor series affords the increment in the 
function: • 

h^ h3 
a 6 



2, According to the definition of the linear 
increment s 



ri = h Y' , n = Ti Y' 
n— 1 n—x ' n— 3 n— s 
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3, Application of the Taylor series results in 
the values of the derivatives in the foregoing intervals 



h'^ h'' IV h* V 

yt = Y« ^ h Y" + ii- Y<" - Y„ + ~ Y . 
n-1 n n 2 » 6' ^ 2.4 



Yi = y«n _ 2hY"„ + -^ ^^'^'"n h*Y^-... 

n-3 n n 2 ^ 6 24 

Yt = Y« - 3h Y" + ~ h®Y«" -~hY„ -+---hY- 

n-3 n ■ ^ 2 ^6 ^24 



4, Suhstitut ion of the previous values in the ex- 
pressions for the linear increments gives 



3 , 4 s 

r, - h Y I - h^Y" + ™ Y I n -. ~ Y^^ + — Y^ - 
^n-i ^^n ^^n+2^n 6 n ^24 n 



■n„ a = h Y' - 2h^ Y" + h Y > " - h Y„ + — h Y, 
''n— s n n 3 n 6 n 24 ^ 



2 . 43 8 4 IV 16 5 J 

■ + h Y>" - ^ h 

3 n 6 

T\ = h Y< - 3h^ Y" + ^ h^Yi" - ^^h'^ Y^^ + -§i h^Y^ 
'n-3 ■ ^ ^ 2 ^ e n 24 ^ 



5. By definition of the first and higher difference 
n— 1 n n— 1 

'n— 2 • 'n— 1 'n— 2 'n 'n— 1 'n— 3 

6, Substitution of, the expressions obtained for 
the linear increments in the differences yields 
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3 IV 



24 



4 IV 36 5 ' 
= h Y „-.iHh Y 



n 



24 



7, Kext it is shown that the true increment in the 
function may "be expressed in ,terms of its linear increments 
and differences: 

For this purpose the coefficients a, P, and Y are com- 
puted "by comparison of the Adams series (107) with the 
Tas'^lor series: 



+ P 



|^h^ht"^ 



' 4 IT 14 5 V 
-hY _ + i^h Y, 
n 24 - ^ 



+ V h ■ Y h y\ 

! 24 ^ 



5 V 



^i! Y^^ + J'- 
24 n 120 ■ n 



This n.av be done hy equating the coefficients in the two 
series of terms of the same degree in h, which then 
yields the following system of eq.uations: 



1 

a, — — 
2 

2 6 
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To d e t e rm i n e t he (3 Q^lf J. c i en.t s^.,t h-e—f i rst-three-e qua t i on s 
are employedj then 

• 2 12 8 



From the fourth equation the ord-er of the computation 
error may be ascertained. In the given case the error 
will be of the order 




By taking, in order to increase the integration accuracy, 
a term with fourth order difference 8A'*riQ_4, the Adams 

equation, after the corresponding computations are carried 
out, takes the following form: 

(108) 

With the aid of the Adams series in the previous form the 
increment in the function (increment of the required in- 
tegral) from point to point can he accurately computed up 
to the terms of the fifth order. Trial computation tables 
have shown that the Adams series in form (108) is suffi- 
ciently accurate ^or integrating the differential equation 
of motion in the boundary layer of the free compressible jet. 

It is seen, from series (l08), that before starting 
the computations it is necessary to have prepared a table 
of values of the function, its first derivative, linear 
increment of the function and firs't;- second, third, and 
fourth differences, respectively, at the fifth, fourth, 
third, second, and first points. In other words, in 
order to start the numerical integration of differential 
equation ( r06') by "the Adam s '"me thod , with t he aid of series 
(108), it is first necessary to compute the integral for 
the above five points. This initial computation is then 
carried out with the aid of the Taylor series on the basis 
of the fact that at the initial point the values of the 
function and its derivatives are known from the boundary 
conditions. 
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5, Integration of th« Differential Equation of the 

■ ? 

Boundary Layer of the Compressible Jet 
The differential equation to "be integrated is 



I'Hi=_p^l+rI'«s + ( 1 06 ) 



with the boundary conditions 

r(0iy = 9i5 2''(cf>i) = l; l*"(cp^)*0; 
PKcPa) = 0; P»(cp2) = 0 



As has been shown in the foregoing the first five 
rows for the numerical -integration table are to be com- 
puted by Taylor series. The computation is made in 
sequence starting vrith point tp^ (the inner limit of 
the boundary layer). Since the argument cp decreases 
from cpj^ to cpg, the integration interval should be 
taken with negative sign. Let: 

Acp = h = - 0. 05 



With subscript o denoting the values of the function 
at F and its' derivatives at F», F", P ' , and F^''" 
at point cpj^ and v/ith subscript n the corresponding 
values at the point cp^ ^ nh , the Taylor series affords 



n 



n+i 



= P„ ' - hF„ " + ^ 



2 ^ 6 ^ 



IV 



> (109^. 



n+i n n 2 ■ " 
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Moreover, the differential -equation gi-ves 



Ft" = r^. , 1 + rl-tg + . a±i (1094) 

n+i n+i n+i- 1+ rF'^ 



and "by further differentiating, 

.IV _ ^■'^ n+1 n+i-'* n+i-" n+i ' n+i^ n+i 



1 +.rr.^^^, • • (1+ rr»s^^.)= 



■pji ji I p II ^ — — 

- ^'n+i ^1 + ^^^'^+1 (1095) 



yr + rj's 

^ n+i 



The . dif f erentiati on ' is limited to the fourth derivative 
in. order to avoid the very complicated expression for 
the fifth derivative. Computations show that the accuracy 
thus obtained is entirely satisfactory. 

The procedure of the preliminary computations is 
as follows. The "boundary conditions give: 

PCcpi) = cp^; i"(q)i) = 1; i"'(tpi) = 0 

With a given value of np^ for given value of the param- 
eter r 

^•■"(q)^) and (cp^) 

are computed. Then hy means of (l07i_3') and (109^_5) 

the values of I", P', P", P ' " , and P^^ at a neighhor— 

ing (second) point (cp^^ — h) are obtained by (107) and 
(109). The same procedure is- followed for all five 
init ial"" points up to (cp^ — 4h). The values of P are 
required only up to point 3. In the remaining rows of 
the table is not required (for the Adams method). 

By denoting the linear increment of the fuixction P and 



I illlllllllli 
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its first two derivatives F', T" by ti = hP'^; 
= hl"'^; = hS'"»jj tables are constructed of 

the linear increments of the functions and their 
differences for the computed five points. The. tables 
of differences are computed by the formulas given in 
section 4: 

In = ^ J"n5 
^ ^n-i -= ^n-™ ^n-i' - 

A Tln-3 ~ ■'in— I + 3t1jj_3 — Tl^—a; 

^\-4 = ^n-4'ln-i-'S^n_3-4^n-=+^n-4 

These formulas are also used for the differences p., 
where, however 

•^n = ^^"n' ^^n = ' "n 
from the Adams series furthermore: . ' 

2 " ^ 12 " ^ 8 720 " 

follow the increments in the function and its derivatives 
at the sixth, seventh, and so forth, points: 

iJ^-'ln^ A^ln-x^--^ "^n-.*! ^' ^n-, ^ ^ 

= '^"^1 ^"-'^fi ^^"-i ^V,*f|i A_ 
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each time increasing the computation table "by one row. 
The ca^lculat ion is "brok-en -of.f.,'when for a given initial 
value of cp^ there is ohtained at the -end of the tahle 
a point cpg satisfying simultaneously the two boundary 
conditions at the outer limit of • the boundary layer: 

PtCcps,) = O;' P"(cpg) = 0 

If one of the derivatives becomes zero ahead of the 
other it indicates that the initial val,ue of. zp^ ..was - 
not we'll chosen, "in this case it is necessary to' 
recompute the table, assuming a new value of cp^^, and 
so forth, until a successful result is obtained. The 
computations carried out show that usually three to 
four, approximations are sufficient to obtain an accurate 
integral table for a given value of Ba. 

In the present paper the computations were carried 
out for air (K = 1,41) first for Ba = 1, Subseq.uently, 
for the other values of Ba, it was possible to obtain 
a good result in a smaller number of trials by interpo- 
lation between the known results for Ba = 1 and Ba ~ 0 
(the latter corresponds to the solution of Tollmien which 
gives the extreme values cp^ - 0.981 and cpg = —2.04), 

By this method of integration the final tables V, 
VI, and VII were obtained for the values of the function 
P and its derivatives F ' , F", and I"'', corresponding 
to various values of the argument cp. Table V was obtained 
for Ba^ = 1.0 (cp^^ = 0.923;: 9^ = - 2.04); table VI for 
Ba^ = 0.9 (cpi = 0.935; cd^ = -2.04); table VII for Bao = 0.5 
(cpi = 0,968; cpg = —2,04). Moreover, for comparison there 
is presented table VIII in which the same magnitudes are 
given for Bao ~ 0 (Tollmien solution), that is, for the 
case of an incompressible fluid. 

Tables V to VIII include also a number of auxiliary 
magnitudes 

_ f; — ; -^') 

by which the longitudinal and transverse velocity com- 
ponents, velocity heads, and so forth, can be computed: 
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t 



P V ■ 

'^o o 



V (pi ' - 



K 



u 



o 

P U 2 
^0 0 



K 



P f2 
K 



6, Pundamental Properties of the Turbulent Boundary Layer 
in the Plane— Parall el Plow of a Compreg sible Gas 

The foregoing solution of the differential equation 
of the boundary. layer of a compressible gas jet makes it 
possible to estimate the effect of the compressibility on 
a turbulent jet of high velocity 

(a) Greometry of the High— Velocity Jet 

As is seen from the tables of numerical integration 
the effect of the compressibility is first of all to 
decrease the thickness of the boundary layer (region of 
mixing of the jet with the surrounding fluid) with in- 
creasing flow velocity; that isj the coordinate (p^ 
of the boundary of the constant velocity core-(u = u^) 
decreases viith increase in Bi., In particular for air 
(K = 1,41) the following values of cpj^ for various 
values of Ba are obtained: 



Ba^ = 1 

Ba 0,9 
o * 

Ba^ =0,5 

Ba 0 (incompress- 
ible fluid) 



cpj^ = .0;923 
cpi- = 0.934 
= 0.966 
(p^ = 0.981 



The relation between cp^ and Ba in terms of the 
factor of compressibility; 



= 2(K - 1) 



Ba^s 



1 + 



K - 1 



Ba, 



(110) 



-NA^A Technical MemOrandiim Ko.- 'l058 71 



can be aijproximat ed by the linear expression' 



cpj^ = O.-9 8'l [1 - 0. 06 r] 



(111) 



In c-QJitrast. to cp^ the ooQ.rcLiiiate of 'the' outer limit cpg 
(u = 0) of the ■'boundary la.y er' .' 'cpg • remains constant with 
change in the flow velocity: 



cpg = constant = — 2,04 



In view of the*' foregoing the nondimenslonal width 
of the boundary layer is connected with the compressi- 
bility factor by th-e relation 



A very important geometrl.c characteristic of the 
jet is the surface separating the initial mass flowing 
under the plane OA, figure 1, frojn the associated mass 
consisting of- the pa.rticles entrained from' the s'lirround— 
ing flow of the gas at .rest. This part itl on ' surf ace 
should also be a flow surface since it, "is the' boundary 
of the initial stream of constant mass floWper second. 
The value of the stream function Si^ ^ the surface 

should be equal to .zeroi 



b = ~- = 3 . 02 [l - 0. 02r ]. ' 



(112) 



But according to eq.uation (,17b), 




(113) 




surface at the ' 



(113b) 



The foregoing expression is thus t.he f.undament'al condition 
for determining the partition surface of the Jet, It is 
seen that the boundary o.f. the core of constant flow is a 
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plane the position of which is determined "by the value 
of the nondinensional coordinate 

, [y] = cp ,(0) (I13c) 

The integration tables in which the values of PC^jr') are 
given afford, by interpolatio=o, the corresponding "raiues: 

Ba. ^ I «p ss - 0, 2 05 

Ba = 0, 9 cp^ = - 0.2 00 

Ba = 0, 5 cp^ « - 0, 192 

Ba 0 ( incompress-' • cp * 0,185 
ible fluid) - 

It is interesting to note that the r-elation "between 

CO and r is also linear.; 
^ 3 

9 = - 0,185 [1 - 0,11 r] (114) 

3 

Moreover,- it has "been showTi that in contrast to the 
■boundary- layer width which decreases with increase in 
the flow velocity, the width of the core of constant 
mass flow increases with inor-ease in velocity. The 
latter effect is explained by the deformation of the 
profile of velocity heads in the boundary layer due to 
the effect of the compressibility of the g&s at high flow 
velocities, 

(b) Velocity Profiles 

Por grea-fcer clearness in comparing the veXocity 
fields obtained for various values of Ba it is neces-sary 
to choose an absolute (independent of Ba) system of co- 
ordinate's. Lay off on the axis of ordinates: 

- CPs 

and on the axis of abscissas 
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Then for all values of Ba the values of the coordinates 
on the inner (cp-^ = 1) and outer (cpg = 0) 'limits of the 
"boundary layer are the same, regardless of the fact that 
the true coordinate of the inner limit cp^ changes with 
change in Ba, • 



Figure 6 gives the curves' F', P ' / for 
along \tfith the curve P' corresponding to the 
ihle fluid (« ss i) for comparison. As may 
velocity profile of the "boundary layer 



Ba = 1 
incompress- 
"be seen the 
hecomes fuller 



with increase in the velocity in the undisturbed region 
( Ba = 1) while the" P» profile becomes less full. In 
general, however, a weak effect of the compressibility 
on the velocity profile of the jet is observed up to 
Ba = 1. 



Figure 9 shows the same curves for Ba = 0.5, The 
character of the curves is the same but the effect of 
the compressibility is now so weak that it may be practi- 
cally neglected. 

Figure 10 gives the density curves in the boundary 
layer for Ba = 0.5 and Ba = 1. 

Figure 11 gives the velocity heads pi^/K for 
Ba = 1 and Ba 0, 

Figure 12 shows a comparison of the transverse 
velocity component distribution: 

■7 _ cp P' - F 



for air at Ba = 1 - and for the incompressible fluid 
(lollmien result), 

(c) Hate of "Flow Discharge 

As established in, a previous part of this paper the 
nondimensional rates of gas flow are subject to formula 
(75) vjhich, applied to the entrained mass, gives 



ms = 



Pa + a 



I- 



Sep 



— dtp 



(115) 
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The nondiinensix)nal 'magnitude of the lni,t ial' mass in the 
jet is - • ' ', 




(116) 



The integral talsles 'give the values i'g(BaQ) and Fj^CBa^) 
and approximate integration ot the values of the integrals 



Ax 



for various values of 
a = 0. C845 it is necessary to determine the magnitudes 
aAg, aA^ and find the values of the entrained and initial 
masses, (Computation shows that the values of aAg and 
aAj^ are relatively very small.) 

It is interesting to note that initial mass, entrained 
mass, and total mass flow through the 'boundary in the 
"boundary layer of the plane-^parallel flow of the high- 
velocity Jet are linear functions of the compressibility 
factor: 

(a) Initial mass: 

= 0, 9 81 [1 ~ 0, 06r ] 

(h) Entrained mass: 

fflg = 0.383 [1 " O^lOr] (lis) 
(c) Total mass rate of flow; 




1 



Bao, Further, assuming with Tollmien 



(117) 
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+ Big = 1,569 [1 - 0. 07 r] - (ll9) 

In the particular cfise of an incompressi"ble • gas the con- 
ventional Tollmien valties' are: 

=,0.981, fflg = 0.388, m^^ + mg 1.369 

(d) Prictional Stresses 

According to ,f ormtila (93) the nondimensicnal f rictional 
stresses at the boundary of the core of -constant mass flow 
is : 




(120) 



With the use of the available integral tables the magnitti.des 
corresponding to the various values of the compressibility 
factor r are obtained by the trapesoidal method. It is 
interesting to note that for the stresses also a linear 
relation, obtains: 

-^-2 = ^3 = 0.27 [1 - 0.167 r] * (121) 



where' 0,27 = is the nondimens ional value of the 

frictional stress in the incompres sible .fluid, 

(e) Conclusions with Regard to the High Velocity Jet 

The laws of variation of th.e fundamental properties 
of the turbulent boundary layer of a plane— parallel stream 
of compressible gas (air) at high flow velocities have 
been obtained in the foregoing. As in the case of the 
"heated" ' jet' the most " interest ing result obtained was 
that the effect of the compressibility on the fundamental 
properties of the high—velocity jet is negligible. This 
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conclusion is valid up to flow velocities attaining 
the velocity of sound. (3a = 1), , 

In particular, on passing from small values of the 
Bairstow numtier ( B.a " 0) to Ba = 1, the angle of di- 
vergence of the boundary lawyer decreases "by '2 percent; 
the angle of dissolution of the core of constant mass 
flow increases hy 11 percent; the nondimensional magni- 
tude of the entrained- mass .de.cr.eases "by 10 percent; the 
nondimensional magnitude of the initial mass decreases 
hy 6 percent; and lastly the nondimensional value of the 
frictional stress at the Ijoundary of the core of constant 
mass flow decreases by 16.7 percent. The effect of the 
c-ompressihility on ■ the properties of the jet is 'so small 
that at Ba = 0,5 it may "be entirely neglected without 
impairment of the results, hence that dynamic similarity 
of the stream prevails over a very wide range of variation 
of the Reynolds num"ber, the temperature a"bove the surround- 
ings and of Ba. The second result in order of imTJortance 
is that the variation of the fundamental properties with 
r is linear, where r is assumed as the c ompres s i'bi li ty 
fact or 



and only the outer "boundary of the Jet remains unchanged, 

V, POSSIBILITY OF DIEBCT APPLICATION OP THE 
PEANDTL LAV/ OF TURBULENT FRICTION TO 
THE CASE OF COMPRESSIBLE JET 

In section I of the present investigation there was 
obtained, for the free Jet of compressible gas, the fric- 
tion law 



which in those cases' where the compressibility- may> be 
neglected (p = constant) reduces to the known Prandtl 
law Obtained fpr incompressible gas: 
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•T = 



(124) 



In comparing eq.uations (123) and (124) t-lie queBtion 
no.turally arises' as, to what extent the accuracy- o-f the 
c ompr es s itle flow investigation is improved "by introducing 
the new friction law and whether this does not lead to 
a considera"ble complication of the problem. Since the 
solution of the prohlem of the boundary layer of a high 
velocity jet was obtained by the new friction law (123,) 
it remains to be explained whether it is possible with- 
out introducing large errors to simplify the given problem 
by direct application of the Prandtl friction law (l24). 
By the latter the differential equation of motion in the case 
of free turbulence is written as follows: 



• bu 5u b 

p u ■ — - + p V — = — 
^ Sx ^ by 5y 



, 2 / Su\ 2 -I 



(125) 



Setting, as pr eviously ^ as sum ed ; 



PKcp) 



pu pv 



= a 



00 '"^00 



(cpP ' - I") ; 1 = cx; 



p ( 



the differential equation is transformed to 



- . = i ^(ir^.^ 
_ .. '■9. .. 



3cp ' 



(126) 



The differential equation in final form for boundary layer 
of a two— dimens ional free jet then is: 
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— I" = p m + 



3 rZii: ^ pnl « I^kJ: (127) 
2 K i K * J K 



As may be -seen, the previous eq.uation obtained by the 
application of the Prandtl friction law is considerably, 
more complicated than eq.uation (20) which was obtained 
on the basis of friction law (123) derived specially for 
free jets of a compressible ga?. In the case of large 
velocities, as shown in section III, a density distri- 
bution is obtained for which 



y~i~+~r?"»2 [l" + yi + rF«2] 
[l+rF'^]^''^ [1 +</! +"r~F '«3 



( 128a) 



(I28b) 



Substitution of the foregoing equations in equation 
(127) leads to the follov/ing form of differential equa- 
tion of the boundary layer of a free high— veloc ity jet: 

1 + ^ y 1 + rP's 

pill-. pyi+rP's+rP's ■ _j. ~-Z-Z ( 129 ) 

[l + rF ' ^ ] [l +y 1 + ri" ' 2 J 



Oomparison of the foregoing equation with the correspond- 
ing equation (l05) gives 

r — ' rPiF"^ 



obtained v/ith the improved friction law (l23) shows con-» 
clusively that the direct extension of the Prandtl fric- 
tion law to the compressible case not only fails to sim- 
plify the study of the latter but on the contrary gives 
a solution considerably more cumbersome than the one ob- 
tained by the more accurate friction law. In view of the 
fact, however, that for a number of aerodynamic problems 
of a compressible gas it may be more convenient to use 
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the ^rei-ctd-o-n --law 4.-n t^he -Pran-dt 1— f ornr, it wa"s* Considered' 
useful for. solving equation ('129). by the Ad^ms method 
and comparing the results with the solution of equation 
(106). 

The following results were ohtained; 

(a) The qu&litajjive ' r esults of the two solutions 
of the problem are the same (with increasing velocity 
the nondimensional values of the frictional- stresses 
decrease, the width of the boundary layer decreases, 
the velocity field is slightly deformed, and the width 
of the core of constant' mass flow increases, 

(b) The CLuant itat i ve results of the two solutions 

differ negligibly. Parti cular'ly -by the use of the 
Prandtl friction law (equation 129) it is found that 
with increase in the Bairstow number from Ba^ 0 to 
Ba^ = 1 .the nondimensional frictional stresses decrease 
by 15 percent; the boundary layer thickness decreases by 
1,3 percent; the divergence angle of the core of constant 
mass flow increases by 16 percent. Correspondingly with 
the ijiprpved friction law (equation (106)) the nondimen- 
sional friction force decreased by 17 percent, the bound- 
ary layer thickness decreased by 2 percent, and the di- 
vergence angle of the constant mass flow core increased 
by 9 percent. 

The application of "the friction law of the incom- 
pressible fluid to the ■ c onpress ible flow thus.lea'ds to 
results that differ slightly from the results .obtained 
with the use of the corrected friction law, In the case 
of the free Jet, however, it is of advantage to apply 
the corrected friction law since it leads to less cumber- 
some, and more readily solvable equations, 

&ENEIIAL COSTCLUSIONS 



■In the present paper ' the ' t he ory of fr ee , turbulence 
and the two— dimensional free Jet was extended to a com— 
.pressihle fluid. In constructing the theory the tur- 
bulence hypothesis of Taylor (vorti.city interchange) was 
used in preference to the Prandtl hypothesis (jnomentum 
interchange). This was done because the former w.as con- 
firmed by test results on the velocity and temperature 
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fields v/hereas the second, while leading to the same 
velocity fields as the first, strongly deviated from 
the- exper iment.ally ohtained temperature fields. 

The differential equation of motion (12) for the 
two— dimensional. ;free .compressible jet leads, in the 
particular case of an incompressi"ble fluid, to the well— * 
known equation of Prandtl-rTollmien-Schlicht ing. 

The boundary layer considered was that of an infi- 
nite^ plane-parallel flow. An infinite flow of this 
kind' assures constancy of the boundary conditions in 
the flow direction; (l) at the outer lijait of the 
boundary layer the flow velocity is equal to zero 
(u = 0)^ (2) at the inner limit the undisturbed flow 
is of constant velocity u = Uq, The constancy of the 
boundary conditions (of the velocities, temperatures, 
densities) at the edges of the boundary layer justifies 
the assumption of aerodynamic similarity (similarity of 
velocity, temperature " fields, etc.) at the various cross- 
sections of the flow, that is, the existence of absolute 
distribution laws of temperatures, velocities, densities, 
frictiona.l stresses, and so forth, as v/as the case with 
incompressible fluids. The absolute differential equa- 
tion o.f motion (20) obtained for the boundary layer of 
an infinite stream is solved for the two cases; 

(1) Jet whose ten-perature differs from that of the 
surroundings and whose velocities are small in comparison 
with the velocity of sound ( 3a ^ 0,5). 

(2) Jet o,f high velocities (up to B'a = 1) the tem- 
perature of which is equal to that of the surroundings. 

The first pase Involved the effect of compressibility 
due to the difference in the temperatures inside and out- 
side the jet; the second case the effect of the compress- 
ibility due to the high flow velocities. 

The fundamental conclusion derived from the results 
of the present investigation is that the effect of the 
compressibility of the fluid on the laws of flov/ in the 
free jet is small. It is also of interest that a lower- 
ing in the temperature of the moving fluid below that of 
the surrounding medium produces fundamentally the same 
effect on the properties of the flow independent of 
whether the lowering was brought about by a cooling of 
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the jet or- by the conv,ersion. of the heat int b- kinet ic 
energy associated with the large fiow^' ve-locit-ieB ; 

In particular: , 

(a) The nondimens ional stresses' oTstained in the air 
jet at moderate velocities and at a teapferature helow 
that of the surroundings hy 6 0° C ivere the same as for 
the jet with large velocities at ; Ba^ =1,0 and have 

a stagnation temperature equal to that of the surround- 
ings, (At Ba =1.0 the local temperature was lower 
than the stagnatioh temperature by about eC.C,') ' 

(b) The nondimensional velocity fields in both* cases 
practically agree (fig. 13), 

(c) The angle of divergence of the. core of constant 
mass flaw in both jets increases by the same amount of 
approximately 11' percent as compared with the "incompress- 
ible jet, 

, Special .compres-sib ility factors are introduced for 
the nonisothermal jet ( S) and for -the 'jet wifh high 
velocities (r) and all fundamental properties of the 
boundary layer were expressed as lineair functions of ' • 
these factors, ' . ' 

The small effect of the c omrir e s s i bi 11 1 y of the fluid 
permits the conclusion that the free jet maintains its 
dynamic similarity over a wide variation of Bairstow 
number and temperature differences as well as for a wide 
range of values of the Reynolds number'; 

The s-uthor has not undertaken the investigation of 
jets escaping from openings of finite diameter and of 
wakes behind bodies in a compressible fluid. "In these 
cases the construction of a theory is effective only at 
great distances from the nozzle (or body) where the vel— ' 
ocity and temperature difference are relatively small 
and hence the effect of . c ompr e s s ib i 1 i ty practically, 
inappreciable though the study of this effect, involves 
great difficulty.. . . • • • _ . ' 

In concluding, the author wishes to express his 
thanks to V, K. Soladkin for his practical a'ssistance 
in the mathematical part of this paper and to A. T, 
Cherniavsky, K. L. Kozlovsky, and L. I. Orlov, who, 
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under the supervision of Solcdkin, carried out all the 
numerical computations, 

V 

Translation Ijy S. Reiss, 
ITational Advisory Committee 
for Aeronautics. 
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TABLE I.- AUXILIiRY mOIIONS OF NONISOTHERII4L JET. 



t 


^0 




F," 


3iF 
S 


3 if' 
S 


3 if" 

s 


0,98 


0,9810 


1.0000 


0 


—0,008 


—0,037 


0,977 


0.90 


0.8988 


0.9977 


0,0746 


-0,002 


-0.118 


0,969 


0,80 


0,8019 


0,9858 


0.1600 


0,019 


-0,213 


0,933 


0,70 


0,7021 


0,9667 


0,2370 


0,044 


—0,304 


0,878 


0.60 


0,6080 


0,9383 


0,3010 


0.079 


-0,391 


0,799 


0,50 


0,5152 


0,9064 


0.3580 


0,122 


-0/166 


0,711 


0,40 


0.4271 


0,8675 


0,4050 


0,173 


—0,531 


0,609 


0,30 


0,3430 


0,8261 


0.4430 


0,228 


-0,584 


0,504 


0,20 


0,2619 


0,7789 


0.4730 


0,287 


-0.629 


0,389 


0,10 


0,1868 


0,7316 


0,4770 


0.365 


-0.634 


0,316 


0 


0,1161 


0.6799 


0,5110 


0.420 


-0,684 


0.16 


-0,10 


0,0508 


0,6297 


0.5190 


0.487 


-ft694 


0,053 


-0,20 


-0.0099 


0,5760 


0,5200 


0.558 


-0,693 


0,001 


—0,30 


-0,0647 


0,5252 


0,5120 


0,627 


-0.682 


-0,090 


-0,40 


-0,1145 


0,4728 


0,4980 


0,693 


-0,662 


-0,174 


-0,50 


-0,1590 


0.4235 


0,4800 


0,757 


-0,633 


—0,246 


-0,60 


-0,1995 


0,3731 


0.4590 


0.818 


-0^95 


-0.313 


-0,70 


-0.2336 


0,3280 


0.4340 


0,924 


-0,505 


-0,402 


-0,80 


-0,2653 


0,2833 


0,4070 


0,924 


-0,505- 


-0,402 


-0,90 


-0,2912 


0,2421 


0,3770 


0,973 


-0.450 


-0,428 


—1,00 


-0,3134 


0,2035 


0,3450 


1,014 


-0,401 


-0,442 


—1.10 


-0,3322 


0,1664 


0,3120 


1.050 


-0,345 


—0,448 


—1,20 


-0,3468 


0,1345 


0,2770 


1,083 


-0.291 


—0,437 


-1,30 


—0,3585 


0,1062 


0.2370 


1,108 


-0338 


—0.405 


-1,40 


-0,3681 


0,0805 


0.2040 


1,132 


-0,196 


-0,376 


-1,50 


—0,3750 


0,0587 


0,1670 


1,144 


—0,140 


-0,329 


-1,60 


-0,3797 


0,0412 


0,1300 


1,155 


-0,102 


-0,269 


-1,70 


-0.38S3 


0,0270 


0.0930 


1,161 


-0,068 


-0.199 


-1,80 


-0,3853 


0,0165 


0,0557 


1,183 


-0,040 


-0,117 


—1.90 


-0,3870 


0.0104 


0,0177 


1,171 


-0.024 


-0,021 


-2.00 


-0,3874 


0,0078 


0,0043 


, 1.178 


-0.016 


0 


-2.04 


-0.3880 


0 


0 


1.181 


— OfiOS 


0 



TABLE II.- JET FUirCTI0ir8 FOR Atg - -60"C. 
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F 


F' 


F" 


X 


fL 


fF'—F 


fF'—F 


F"- 

X 


e 


X 


1,005 


1,005 


1,0 


0 


1.0 


1,0 


0 


0 


1,0 


1 


0,9 


0,8989 


0,9951 


0,0985 


0,9939 


1,0012 


—0,0040 


-0,0040 


0,9963 


0,9754 


0,8 


0,8025 


0,9807 


0,1836 


0.9866 


1.9941 


-0,0179 


-0,0181 


0,9750 


0,9458 


0,7 


0,7031 


0,9595 


0,2586 


0,9793 


0,9798 


-0,0315 


—0,0322 


0,9401 


0,9153 


0,& 


0,6099 


0,9294 


0,3207 


0,9721 


0,9560 


—0,0523 


-0,0538 


0,8885 


0,8853 


0.5 


0,5180 


0,8945 


0.3755 


0,9650 


0.9270 


-0.0708 


—0,0734 


0,8291 


0,8549 


0,4 


0,4313 


0,8549 


0,3200 


0,9579 


0,8925 


-0,0893 


-0,0932 


0,7630 


0,8242 


0,3 


0,3486 


0,8166 


0,4554 


0,9509 


0,8535 


-0.1051 


- 0,110s 


0,6927 


0,7934 


0,2 


0,2691 


0,7635 


0.4826 


0,9439 


0,8089 


-0,1164 


-0,1233 


0,6176 


0,7622 


0,1 


0.1960 


0,7164 


0,4848 


0,9369 


0.7646 


-0,1244 


—0,1328 


0,5478 


0,7305 


0 


0,1264 


0,6631 


0,5149 


0,93 


0,7130 


-0,1264 


—0,1359 


0,4728 


0,6989 


-0,1 


0,0629 


0,6129 


0,5203 


0,9232 


0,6639 


—0,1242 


-0,1345 


0,4068 


0,6674 


-0,2 


0,0039 


0,5589 


0,5200 


0,9164 


0,6099 


-0.U55 


-0,1262 


0,3410 


0,6351 


-0,3 


—0,0496 


0,5082 


0,5098 


0,9096 


0,5587 


-0.1028 


-0,1130 


0,2840 


0,6026 


-0,4 


-0,0979 


0,4567 


0,4937 


0.9030 


0.5058 


-0,0848 


-0,0039 


0,2310 


0,5704 


-0,5 


—0,1403 


0,4084 


0,4739 


0,8963 


0,4557 


-0,0639 


-0,0713 


0,1861 


0,5371 


-0,6 


-0,1788 


0,3583 


0,4513 


0,8898 


0,4028 


—0.0362 


-0,0407 


0.1443 


0,5042 


-0,7 


-0,2125 


0,3143 


0.4251 


0,8832 


0,3559 


-0,0075 


-0,0085 


0,1108 


0,4711 


-0,8 


—0,2422 


0,2705 


0,3971 


0.8767 


0.3086 


0,0258 


0,0294 


0,0835 


0,4375 


-0,9 


-0,2670 


0,2308 


0,3664 


0,8702 


0.2653 


0,0593 


0,0681 


0.0612 


0,4035 


-1,0 


—0,2880 


0,1941 


0,3341 


0,8639 


0.2247 


0,0939 


0,1087 


0,0436 


0,3698 


-1,1 


-0,3061 


0,1575 


0,3009 


0,8575 


0,1837 


0,1329 


0,1550 


0,0289 


0.3353 


—1,2 


-0,3203 


0,1268 


0,2662 


0,8512 


0,1490 


0,1681 


0,1975 


0,0189 


0,3008 


-1,3 


-0,3317 


0,1001 


0,2270 


0,8449 


0.U85 


0.2015 


0,2385 


0.0118 


0,2657 


-1,4 


-0,3401 


0,0757 


0.1947 


0,8387 


0.0902 


0.2349 


0.2801 


0.0068 


0,2308 


-1,5 


—0,3468 


0,0552 


0,1589 


0,8326 


0,0664 


0,2639 


0,3170 


0,0037 


0,1958 


-1,6 


-0,3515 


0,0387 


0,1234 


0,8264 


0,0468 


0,2896 


0,3504 


0,0018 


0,1597 


-1,7 


—0,3544 


0,0254 


0,0881 


0,8203 


0.0310 


0,3112 


0,3794 


0,0008 


0.1237 


-1,8 


—0.3558 


0,0155 


0.0528 


0,8143 


0,0190 


0.3279 


■ 0,4027 


0,0003 


0.0878 


-1,9 


-0,3581 


0,0081 


0,0172 


0,8083 


0,0100 


0,3427 


0,4249 


0,0001 


0,0513 


—2,0 


-0.3581 


0,0020 


0 


0,8024 


0,0025 


0,3541 


0,4414 


0,0001 


0,0150 


'-2,04 


-0,3581 


0 


0 


0,8000 


0 


0,3581 


0,4475 


0 


0 
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TABLE III.- JET FOITCTIONS FOR & t « 60°C. 
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0,96 


0,96 


1,0 


0 


1,0 


1.0 


0 


0 


1.0 


1,0 


0,9 


0,8990 


1,0 


0,0537 


1.0037 


0,9963 


0.001 


0,0010 


0,9964 


0.9781 


0,8 


0,8016 


0,9903 


0,1398 


1,0098 


0.9807 


-0,0094 


-0,0093 


0,9712 


0,9416 


0,7 


0,7011 


0,9731 


0,2180 


1.0159 


0.9579 


- 0,0199 


-0,0196 


0.9321 


0.9060 


0,6 


0,6064 


0,9469 


0,2837 


1,0221 


0,9264 


-0,0383 


-0.375 


0,8772 


0,8703 


0,5 


0,5125 


0,9154 


0,3426 


1,0284 


0.8901 


-0,0548 


-0,0533 


0,8149 


0.8342 


0,4 


0,4235 


0,8787 


0,3918 


1,0340 


0,8498 


—0,0720 


—0,0696 


0,7467 


0,8027 


0,3 


0,3384 


0.8378 


0,432 


1,0409 


0,8049 


-0,0871 


—0,0837 


0.6743 


0,7642 


0,2 


0,2562 


0,7917 


0,4646 


1,0500 


0,7540 


—0,0979 


-0,0932 


0.5970 


0,7143 


0,1 


0,1796 


0,7448 


0,4702 


1,0537 


0,7068 


-0,1051 


—0.0997 


0.5264 


0,6944 


0 


0,1075 


0,6938 


0,5075 


1,0600 


0,6545 


-0,1075 


-0,1014 


0,4542 


0,6604 


-0,1 


0,0411 
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0,6041 


-0,1055 
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-0,0968 
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- 0,3 
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0,5139 
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-0,0840 
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0,5577 


-0,4 
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0,4864 


0,5018 


1,0861 


0,4478 
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-0,0604 


0,2178 


0,5244 


1-0,5 


-0,1743 


0,4368 


0,4853 


1.0928 


0,3997 


—0,0441 


—0.0404 


0,1746 


0.4904 


1 

1-0,6 


-0,2155 


0,3850 


0,4662 


1,0994 


0,3502 


—0.0155 


—0.0141 


0,1348 


0,4573 


- 0,7 


-0,2516 


0,3392 


0,4418 


1,1061 


0,3067 


0.0142 


+0,0128 


0,1041 


0,4245 


0,8 


-0,2836 


0,2932 


0,4157 


1.1130 


0,2634 


0,0490 


0,0440 


0.0773 


0,3908 


-0,9 


-0,3106 


0,2512 


0,3863 


1,1197 


0,2244 


0,0845 


0,0755 


0,0564 


0,3586 


-1,0 


-0,3334 


0,2121 


0,354b 


1,1265 


0,1883 


0,1213 


0,1077 


0,0399 


0,3262 


-1,1 


-0.3532 


0,1730 


03217 


1.1334 


0,1526 


0,1629 


0.1437 


0.0264 


0,2937 


— 1,2 


-0,3688 


0,1399 


0,2676 


1,1403 


0,1227 


0.2009 


0,1762 


0.0172 


0,2617 


-1,3-^ 


' 0,38)3 


0,1108 ' 


0,2283 


1.1472 


0.0966 


0,2373 


0,2069 


0.0107 


0,2302 


-1,4 


0,3908 


0,0844 


0,1959 


1,1542 


0,0732 


0,2726 


0,2362 


0,0062 


0,1984 


-1,5 


0,3981 


0,0615 


0,1741 


1,1613 


0,0530 


0,3058 


0,2633 


0,0033 


0,1666 


-1,6 


0,4033 


0,0433 


0,1358 


1.1683 


0,0370 


0,3341 


0,2860 


0,0016 


0,1356 


-1,7 


0,4064 


0,0285 


0,0973 


1.1755 


0,0242 


0,3580 


0,3045 


0,0007 


0,1042 


-1,8 


0,4089 


0,0173 


0,0582 


1.1826 


0,0146 


0,3776 


0,3193 


0,0003 


0,0736 


-1,9 


0,4106 


0,0109 


0,0182 


1,1892 


0,0092 


0,3899 


0,3279 


0,0001 
• 


0,0449 


- -2,0 


0,4106 


0,0081 


0 


1,1971 


0,0068 


0,3943 


0,3294 


0,0001 


0,0121 


-2,04 


0,4106 


0 


0 


1,2000 


0 


0.4106 


0,3420 


0 


0 



TABLE IV.- JET FUNCTIONS FOR ttQ = 120°C 



0,942 
0,9 
0,8 
0,7 
0,6 
0,5 
0,4 
0,3 
0,2 
0,1 
0 

-0,1 
-0,2 
-0,3 
-0,4 
-0,5 
-0,6 
-0,7 
-0,8 
-0,9 
-1,0 
-1,1 
-1,2 
-1,3 
-1,4 
-1,5 
-1,6 
-1,7 
-1,8 
—1,9 
-2,0 
— 2,04 



0,942 
0,8991 
0,8013 
0,7004 
0,6051 
0,5105 
0,4206 
0,3345 
0,2513 
0,1734 
0,1004 
0,0328 
-0,0306 
-0,0883 
-0,1408 
-0,1871 
—0,2294 
-0,2664 
-0,2993 
-0.3272 
—0,3507 
—0,3710 
-0,3872 
-0,4002 
-0,4101 
-0,4175 
—0,4229 
-0,4261 
—0,4290 
—0,4305 
-0,4305 
-0,4305 



1,0 
1.0 

0,9939 
0,9783 
0,9535 
0,9233 
0,8877 
0,8477 
0,8024 
0,7556 
0,7054 
0,6558 
0,6018 
0,5503 
0,4976 
0,4475 
0,3951 
0,3486 
0,3018 
0,2589 
0,2189 
0,1788 
0,1448 
0,1148 
0,0878 
0,0639 
0,0450 
0,0296 
0,0180 
0,0113 
0,0084 
0 



0 

0,0386 
0,1253 
0,2044 
0,2713 
0,3316 
0,3824 
0,4243 
0,4586 
0,4652 
0,5051 
0,5170 
0,5200 
0,5153 
0,5045 
0,4891 
0,4707 
0,4474 
0,4220 
0,3929 
0,3614 
0,3287 
0,2932 
0,2520 
0,2180 
0,1792 
0,1400 
0,1004 
0,0601 
0,0185 
0 

0 . 



1,0 

1,0048 

1.0162 

1,0272 

1,0393 

1,0511 

1,0630 

1,0751 

1,0873 

1,0996 

1,1121 

1,1248 

1,1375 

1.1504 

1,1635 

1,1767 

1,1900 

1,2035 

1,2172 

1,2310 

1,2450 

1,2600 

1,2740 

1,2880 

1,3026 

1,3320 

1,3474 

1,3627 

1,3782 

1,3938 

1,4000 

1,4000 



1,0 

0,9952 
0,9781 
0,9524 
0,9174 

0,8784 
0,8351 
0,7885 
0.7380 
0,6872 
0,6343 
0,5830 
0,5291 
0,4784 
0,4277 
0,3803 
0,3320 
0.2897 
0,2480 
0.2103 
0.1758 
0.1419 
0,1137 
0,0891 
0,0674 
0,0480 
0,0334 
0,0217 
0,0131 
0.0081 
0,0060 
0 



fF' — F 

0 

0.0009 
—0,0062 
-0,0156 
—0,0330 
-0,0489 
-0,0655 
-0,0802 
—0,0908 
—0,0978 
—0,1004 
-0,0984 
-0.0898 
—0,0767 
-0,0582 
-0,0367 
—0.0077 
0,0224 
0,0579 
0,0942 
0,1318 
0,1743 
0,2134 
0,25J0 
0,2872 
0,3217 
0,3509 
0.3758 
0.3966 
0,4090 
0,4137 
0,4305 



I 

0 

0,0009 
—0,0061 
—0,0152 
-0,0318 
-0,0465 
-0,0616 
—0,0746 
—0,0835 
—0,0889 
-0,0903 
-0,0375 
-0,0789 
-0,0667 
-0,0500 
-0,0312 
0,0065 
0,0186 
0,0476 
0,0765 
0,1059 
0,1383 
0,1675 
0,1949 
0,2205 
0,2415 
0,2604 
0,2758 
0,2878 
0,2935 
0,2955 
0,3080 



7, 

1,0 

0,9952 
0,9721 
0,9318 
0,8748 
0,8111 
0.7413 
0,6684 
0,5921 
0,5192 
0,4474 
0,3824 
0,3140 
0,2632 
0,2128 
0,1702 
0,1312 
0,1010 
0,0748 
0,0544 
0,0385 
0,0254 
0,0165 
0,0102 
0,0059 
0,0031 
0,0015 
0,0007 
0,0002 
0,0001 
0,0001 
0 



jTABU: t.- BASIC FUNCTZ01I8 FOR Bao - 1.0. 



[i/n 


9 


F 








0 


0,923 


0,923000 


1.000000 


0 


—1,292200 


1 


0,873 


0,873027 


0,998414 


0.062860 


—1,219354 


2 


0,823 


0,823210 


0,993778 


0,121941 


-1,141705 


3 


0,773 


0,773697 


0.986378 


0,177036 


—1,060580 


4 


0,723 


0,724626 


0,976143 


0,228006 


-0,977252 


5 


0,673 


0,676124 


0,963557 


0,274762 


—0,892915 


6 


0,623 


0,628308 


0,948737 


0,317299- 


-0.808650 


7 


0,573 


0,581285 


0,93189& 


0,355643 


-0,725433 


8 


0,523 


0,535149 


0,913241 


0.389875 


—0,644099 


9 


0,473 


0,489987 


0,892973 


0,420096 


-0,565365 


10 


0,423 


0,445873 


0,871294 


0,446463 


—0.489797 


11 


0,373 


0,402878 


0,848388 


0,469138 


—0,417885 


12 


0,323 


0,361055 


0,824438 


0,488317 


—0,349931 


13 


0,273 


0,320447 


0,799612 


0,504203 


—0,286158 


14 


0,223 


0,281103 


0,774070 


0,517004 


—0,226710 


15 


0,173 


0,243049 


0,747960 


0,526947 


-0,171613 


16 


0,123 


0,206313 


0,721421 


0.934237 


-0,120842 


17 


0,073 


0,170912 


0,694578 


0,539099 


-0,074289 


18 


0,023 


0,136859 


0,667548 


0,541734 


-0,031810 


19 


—0,027 


0,104160 


0.640437 


0,542346 


+0,006783 


20 


—0,077 


0,072813 


0,613346 


0,541115 


0,041721 


21 


—0,127 


0,042824 


0,586354 


0,538226 


0,073216 


22 


—0,177 


0,014173 


0,559549 


0,533845 


0,101545 


23 


—0,227 


—0,013136 


0.532993 


0.528120 


0,126950 


24 


—0,277 


—0,039131 


0,506756 


0.521191 


0,149703 


25 


—0,327 


—0,063817 


0,480891 


0.513191 


0,170042 


26 


—0,377 


—0,087227 


0/455451 


0.504225 


0.188220 


27 


—0,427 


—0,109373 


0,430182 


0,494401 


0,204462 


28 


—0,477 


—0,130284 


0,406020 


0,483808 


0,218985 


29 


—0,527 


-0,149985 


0,382111 


0,472528 


0,231989 


30 


—0,577 


—0,168506 


0,358778 


0,460632 


0,243657 


31' 


-0,677 


-0,202122 


0,313971 


0,435229 


0,263629 


32 


—0,777 


—0,231386 


0,271 793. 


0,408024 


0,280012 


33 


-0,877 


-0,256572 


0,232417 


0,379320 


0,293660 


34 


-0.987 


—0,277970 


0,195973 


0,349363 


0,305194 


35 


—1,077 


-0,295872 


0,162582 


0,318342 


0,315029 


36 


—1,177 


-0,310593 


0,132336 


0.286407 


0,323442 


37 


—1,277 


-0,322447 


0,105324 


0.253699 


0,330596 


38 


-1,377 


—0,331764 


0.081618 


0.220329 


0,336603 


39 


— 1,477 


—0,338884 


0,061276 


0,186417 


0.341531 


40 


-1,577 


-0,344135 


0.044348 


0,152059 


0,345442 


41 


—1,677 


—0,347868 


0,030877 


0,117360 


0,348435 


42 


-1,777 


-0,350428 


0,020885 


0,082399 


0,350637 


43 


—1,877 


-0,352159 


0.014402 


0,047254 


0.352225 


44 


-1,977 


-0,353436 


0.011440 


0,01 1970 




45 


-2,04 


-O.35400O 


0 


0 





TABLI .V( a). - AUXILIARY FtMCriONS FOR Bao " 1«0« 



n/n 


T 


It 


IL 

X 


i 


fF--F 

X 


ft 

X 


0 


0.923 


1.00 


1,00 


0 


0 


1,000000 


1 


0.823 


0.9982 


0,9956 


-0,005331 


—0,005341 


0,985817 


2 


0,723 


0,9932 


0.9888 


-0,018875 


-0,019004 


0.946376 


3 


0,623 


0,9856 


0.9626 


—0,037245 


—0,037789 


0,887141 


4 


0,523 


0,9758 


0,9358 


—0,057524 


-0,058951 


0,913826 


5 


0,423 


0,9645 


0.9034 


-0,077316. 


—0.080162 


0,732153 


6 


0,323 


0,9523 


0.8657 


-0,094762 


-0,099509 


0,647276 


7 


0,223 


0.9396 


0.8239 


-0,108485 


—0,115459 


0,562993 


8 


0,123 


0,9269 


0,7783 


-0,117578 


- 0,126851 


0,482403 


9 


0,023 


0.9145 


0,7299 


—0.121505 


-0,132865 


0,407519 


10 


—0,077 


0.9028 


0,6793 


—0,120041 


-0,132965 


0,339627 


11 


-0,177 


0,8918 


0,6274 


—0,113213 


—0.126949 


0,279218 


12 


— OJ277 


0,8819 


0,5747 


—0,101240 


-0,114796 


0,226474 


13 


-0,377 


0,8729 


0,5217 


-0,084478 


-0,096779 


0,181071 


14 


-0,477 


0,8651 


0,4693 


—0,063388 


-0,073272 


0.142613 


15 


-0.577 


0,8567 


0,4188 


-0,038509 


-0,044950 


0,1 10276 


16 


-0,677 


0,8526 


0,3682 


-0,010436 


-0.012240 


0,084048 


17 


—0,777 


0,8479 


0,3206 


-fO.020203 


-1-0,023837 


0.062635 


18 


—0.877 


0,8440 


0,2754 


0,052742 


0,062491 


0,045591 


19 


-0.977 


0,8409 


0,2331 


0,086504 


0,102871 


0,032295 


20 


-1,077 


0,8386 


0,1939 


0,120771 


0,144015 


0.032167 


21 


-1,177 


0,8368 


0,1581 


0,154834 


0,185031 


0,014655 


22 


-1.277 


0,8355 


0,1260 


0.187948 


0,224953 


0,009268 


23 


—1,377 


0,8347 


0,09778 


0.219376 


0,252820 


0,005560 


24 


—1,477 


0,8341 


0,07346 


0.248379 


0,297781 


0,003132 


25 


— 1,577 


0,8337 


0,05320 


0,274195 


0,328893 


0,001640 


26 


-1,677 


0,8334 


0,03706 


0.296087 


0.355276 


0,000794 


27 


— 1,777 


0,8334 


0,02506 


0,313315 


0,375948 


0,000363 


28 


-1,877 


0,8334 


0,01728 


0.325126 


0,390134 


0.000172 


29 


-1,977 


0,8340 


0,01370 


0,326000 


0,392000 


0.000160 


30 


-2,04 


0,8340 


0 


0,326009 


0,392000 


° i 



00 



TABLI VI.- BASIC FUNOTIOHS FOR Bag = 0.9. 



TABLE VI(a).- AUXILIARY FUNCTIONS FOR BBq ° 0.9. 



n,ii 




F 


F' 






0 


0,934 


0,934000 


1.00 


0 


— 1,236624 


1 


0,884 


0,884026 


0,998482 


0,060176 


— 1,168140 


2 


0,834 


0,834201 


0.994042 


0,116818 


-1,095812 


3 


0,784 


0,784668 


0.986862 


0.169762 


—1,02071 1 


4 


0,734 


0,73S558 


0,977130 


0,218895 


- 0,943871 


5 


0,684 


0,686996 


0,965037 


0,264152 


-0,866263 


6 


0,634 


0,639090 


0,950780 


0.305526 


-0,788765 


7 


0,584 


0,591950 


0,934550 


0343044 


-0,712178 


8 


0,534 


0,545664 


0,916538 


0,376771 


-0,637192 


<J 


0.484 


0,500322 


0.896933 


0,406801 


-0.564401 


10 


0,434 


0,455993 


0,875919 


0,433255 


—0,494290 


11 


0,384 


0,412751 


0,853663 


0,456280 


—0,427252 


12 


0,334 


0,370647 


0.830344 


0,476038 


-0,363573 


13 


0,284 


0,329732 


0,806112 


0,492695 


-0,303458 


14 


0,234 


0,290048 


0,781122 


0.506443 


-0,247031 


15 


0,184 


0,251632 


0,729418 


0,525938 


-0.145390 


IB 


0,134 


0.214S04 


0,729418 


0,525938 


-0,145390 


17 


0,084 


0.178695 


0.702959 


0,532063 


-0,100106 


18 


0,034 


0,144212 


0,676249 


0.536007 


—0,058384 


19 


-0,016 


0,1 10720 


0,549393 


0,537957 


—0,020086 


20 


-0,066 


0.079271 


0,622483 


0,538069 


+0,014958 


21 


-0,116 


0,0-18823 


0,595616 


0,536510 


0.046922 


22 


-0,166 


0,019708 


0,568866 


0,533426 


0.07602G 


23 


-0,216 


-0,008067 


0,542296 


0,528951 


0,102454 


24 


-0,266 


-0,034525 


0.515982 


0,523218 


0,126430 


25 


-0,316 


-0,059671 


0,489992 


0,516343 


0,148149 


26 


-0,366 


-0,083531 


0,464367 


0,508438 


0,167823 


27 


-0,416 


—0,106113 


0,439162 


a499591 


0.185627 


28 


—0,466 


-0,127454 


0,414422 


0,489903 


0,201764 


29 


-0,516 


—0,147566 


0.3901&5 


0,479441 


0,216385 


30 


-0,566 


—0,166480 


0,366489 


0,468286 


0,229655 


31 


-0,666 


-0,200830 


0,320846 


0,444128 


0,252695 


32 


-0,766 


—0,230735 


0,277733 


0,417872 


0,271854 


33 


-0,866 


-0,256465 


0,237331 


0,389862 


0,287904 


34 


—0,966 


-0,278298 


0,199809 


0,360375 


0,301420 


35 


-1,066 


-0,296528 


0,165298 


0,329649 


0,312816 


36 


-1,166 


-0,311461 


0,133915 


0,297876 


0,322385 


37 


-1,266 


-0,323418 


0,105750 


0.265232 


0,330332 


38 


—1,366 


—0,332721 


0,080892 


0.231861 


0,336799 


39 


-1,466 


—0,339710 


0.059398 


0,197918 


0,341908 


40 


—1.566 


-0,344717 


0,041321 


0,163524 


0,345769 


41 


-1,666 


-0,348088 


0,026705 


0.128802 


0,348505 


42 


—1,76b 


-0,350173 


0,015568 


0,093859 


0,350308 


43 


—1,866 


—0,351319 


0,007935 


0,058765 


0.351348 


44 


—1,966 


—0,351876 


0.003220 


0,023606 


0,351880 


45 


—2,04 


-0,352000 


0 


0 







? 


X 


F' 


<pF' —F 


fF' — F 




n'n 








X 


X 


0 


0,934 


1,00 


1,00 


0 


0 


1.0000 


1 


0,834 


09985 


0,9955 


—0,005170 


—0,005178 


0,9896 


2 


0,734 


0,9944 


0,9826 


-0.018345 


-0,018448 


0,9602 


3 


0.634 


0,9381 


0,9622 


-a03629S 


-0.036732 


0,9149 


4 


0,534 


0,9801 


0,9351 


-0,056233 


—0,057375 


0,8571 


5 


0,434 


0,9708 


0,9023 


-0,075844 


-0,078125 


07903 


6 


0,334 


0,9606 


0,8644 


-0,093312 


-0,097139 


0,7178 


7 


0,234 


0,9501 


0.8221 


-0,107265 


— aU2899 


a6422 


8 


0,134 


0.9395 


0,7764 


-0,116762 


—0,124281 


0,5663 


9 


0,034 


0.9292 


0,7278 


-0,121220 


-0,130456 


0,4922 


10 


—0,066 


0,9193 


0,6771 


-0,120355 


-0,130920 


0,4215 


11 


-0,166 


0,9102 


0,6250 


—0,114138 


—0,125399 


0,3555 


12 


-0,266 


0,9017 


0,5722 


—0,102726 


—0,113925 


0,2953 


13 


-0,366 


0.8942 


0,5193 


— 0,086427 


-0,096653 


0,2412 


14 


—0,466 


0,8876 


0,4669 


-0,065667 


- 0,073983 


0,1935 


15 


-0,566 


0,8818 


0,41S6 


-0,040953 


-0,046443 


0,1523 


16 


-0,666 


0,8770 


0.3658 


—0,012853 


—0,014656 


0,1174 


17 


-0,766 


0,8789 


0,3182 


40,017992 


+0,020621 


0,08837 


18 


-0,866 


0,8696 


0,2729 


0,050936 


0,058574 


0,06477 


19 


—0,966 


0,8670 


0.2305 


0,085283 


0,098366 


0,04605 


20 


— 1.066 


0,8650 


0,1911 


0,120320 


0,139099 


0,03159 


21 


-1,166 


0,8635 


0,1551 


0,155316 


0,179868 


0,02077 


22 


-1,266 


0,8624 


0,1226 


0,189539 


0,219781 


0,01297 


23 


— 1,366 


0,8616 


0,09389 


0.222223 


0.257919 


0,00759 


24 


—1.466 


0,8612 


0,06897 


0,252633 


0,293350 


0.00409 


25 


-1,566 


0,8609 


0,04800 


0,280008 


0,325250 


0,00198 


26 


— 1,566 


0,8607 


0.03103 


0,303598 


0,352734 


0,00082 


27 


— 1,766 


0,8606 


0,01809 


0,322680 


0,374948 


0.00028 


28 


-1,866 


0.8606 


0,00922 


0,336512 


0,391020 


0,00007 


29 


—1,966 


0,8606 


0.00443 


0,346000 


0,402000 


0 


30 


-2,040 


0,8606 


0 


0.352000 


0,409000 


0 



TIBLI VII.- BASIC FiraCTIOHS FOR Bsq = 0.5. 



TABLE Vll(a).- AUXILIARY JUMCTIQNS FOR Ba^ -0.5. 



n/n 


f 


F 


F- 




P- 


0 


0,966 


0,966 


1,0 


0 


—1,062600 


1 


0,916 


0,916032 


0,998695 


0,051755 


—1,006931 


2 


0,866 


0,866173 


0,994871 


0,100695 


-fl,950190 


3 


0,816 


0,816575 


0,988672 


0,146776 


-0,892768 


4 


0,766 


0,767343 


0,980241 


0,189973 


—0,835025 


5 


0,716 


0,718586 


0,969722 


0,230282 


—0,777326 


6 


0,666 


0,670403 


0,957260 


0,267710 


—0,719989 


7 


0,616 


0,622889 


0,942995 


0,302289 


—0,66331 1 


8 


0,566 


0,576129 


0,927077 


0,334059 


-0,607562 


9 


0,516 


0,530205 


0,909634 


0,363064 


—0,552974 


10 


0,466 


0,485188 


0,890816 


0,389380 


-0,409746 


] 1 


0,416 


0,441143 


0,870741 


0,413065 


—0,448052 


12 


0,366 


0,398132 


0,849551 


0,434213 


—0,398029 


13 


0,316 


0,356202 


0,827360 


0,452897 


—0,349786 




0,266 


0,315409 


0,804300 


0,469221 


—0,303412 


15 


0,216 


0,275784 


0,780477 


0,48327? 


—0,258955 


16 


0,166 


0,237371 


0,756007 


0,495148 


—0,216460 


17 


0,1 16 


0,200195 


0,731000 


0^04952 


—0,175936 


18 


O.Ofab 


0,154277 


0,705543 


0,512776 


—0,137379 


19 


0.016 


0,129646 


0,6?9751 


0,512776 


-0,137379 


20 


—0,034 


0,096306 


0,653702 


0,522886 


—0,066092 


21 


—0,084 


0,064280 


0,627489 


0,525362 


—0,033293 


22 


-0,134 


0,033562 


0.601192 


0,526245 


—0,002316 


23 


— 0,184 


0,004158 


0,574889 


0,525623 


0,026889 


24 


—0,234 


0,023926 


0,548653 


0,523585 


0,054379 


25 


—0,284 


0,050710 


0,522553 


0,520210 


0,080230 


26 


— 0,334 


0,076186 


0,496554 


0,515585 


0,104494 


21 


— 0,384 


0,100378 


0,471015 


0,509786 


0,127252 


23 


— 0,434 


0,123294 


0,445696 


0,502882 


0,148560 


29 


— 0,484 


0,M4954 


0,420745 


0,494952 


0,168493 


30 


— 0,534 


0,165375 


0,396214 


0,486054 


0,187109 


31 


— 0,634 


0,202595 


0,348603 


0,465627 


0,220641 


32 


— 0,734 


0,235169 


0,303191 


0,442078 


0,249635 


33 


-0,834 


0,263319 


0,260276 


0,415834 


0,274503 


34 


— 0,934 


0,287316 


0,220102 


0,387301 


0,295637 


35 


— 1,034 


0,307438 


0,182877 


0,356818 


0,313372 




— 1,134 


0,323997 


0,148791 


0,324729 


0,328029 


37 


— 1,234 


0,337306 


0,1 17976 


0.291307 


0,339895 


38 


— 1,334 


0,347705 


0,090562 


0,256831 


0,349256 


39 


— 1,434 


0,355534 


0,066637 


0.221531 


0,356386 


40 


—1.534 


0,361152 


0,046275 


0,185618 


0,361568 


41 




U,oo717o 


0,016430 


0,1 12640 


0,367232 


42 


-1.834 


0,368320 


0,007004 


0.075854 


0,368330 


43 


-1,934 


0,368702 


0,001261 


0.038998 


0,368702 


44 


-2,04 


0,369000 


0 


0 





n n 


T 




£1 

X 


tF' -F 


fF' — F 
n 


£2 

X 


0 


0,966 


1,00 


1.00 


0 


0 


1,0000 


1 


0,866 


0,9996 


0,0954 


-0,004573. 


-0,004575 


0.9902 


2 


0,766 


0,9983 


0,9819 


-«,0I6S43 




0,9625 


3 


0,666 


0,9963 


0,9609 


-0,032803 


-0.032925 


0.9197 


4 


0,566 


0,9936 


0,9331 


—0,051429 


-0.051760 


0,8650 


5 


0,466 


0,9903 


0,8993 


-0,070088 


- 0.070775 


0.3013 


6 


0,366 


0,9872 


0.8607 


-0,087132 


-a088262 


0,7311 


7 


0,266 


0,9838 


0,8176 


-0,101509 


-0.103181 


0,6575 


8 


0,166 


0,9802 


0,7713 


-0,111871 


—0,114131 


0.5831 


9 


0,066 


0,9767 


0,7223 


—0,117677 


-0,120484 


0,5097 


10 


—0,034 


0,9733 


0,6717 


-0,118506 


-0,121757 


0,4390 


11 


-0,134 


0,9702 


0,6197 


-0.114162 


—0,1 17668 


0,3735 


12 


-0,234 


0,9672 


0,5672 


-0,104474 


-0,108017 


0.3112 


13 


- 0,334 


0,9642 


0,5148 


-O.08971* 


-0.093006 


0.2558 


14 


-0,434 


0,9622 


0,4632 


-0,070106 


-0.072860 


0,2064 


15 


-0,534 


0,9603 


0,4127 


-0,046225 


-0,048136 


0,1635 


16 


-0,634 


0,9584 


0,3637 


—0,013405 


—0,019204 


0,1268 


17 


-0,734 


0,9570 


0,3168 


-'-0,012669 


-i-0,013238 


0,0960 


18 


—0,834 


0,9557 


0,2724 


0,046219 


0,048361 


0,0708 


19 


-0,934 


0,9549 


0,2305 


0,081716 


0,085575 


0,0507 


20 


-1,034 


0,9541 


0,1917 


0,118338 


0,124031 


0,0350 


21 


- 1,134 


0,9535 


0,1552 


0,155297 


0,162870 


0,0232 


22 


- 1,234 


0,9532 


0,1238 


0.191706 


0,201118 


0,0146 


23 


— 1,334 


0,9528 


0,0951 


0,226805 


0,238041 


0.0086 


24 


— 1,434 


0,9527 


0,0699 


0,260034 


0,272944 


0,0046 


25 


-1,534 


0,9525 


0,0487 


0.290152 


0,304621 


0.0022 


2fa 


— 1,634 


0,9525 


0,03098 


0,316710 


0,332504 


0,0009 


27 


— 1,734 


0,9524 


0,01720 


0,338778 


0,355710 


0,0002 


28 


-1,834 


0,9524 


0,00740 


0,355520 


0,373289 


0,0000 


29 


-1,934 


0,9524 


0,00140 


0,366202 


0,384504 


0,0000 


30 


-2,040 


0,9524 


0 


0.3690 


0,38750 j 


0 
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O 

B" 
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P 
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00 
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TABLE VIII 
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BASIC FUNCTIONS FOR INCOMPRESSIBLE JET. 



^6 
n,n 


•P 1 P 


F' 


F" 


/•"' 



0 


0,98 


1 


0,93 


2 


0,88 


3 


0,83 


4 


0.78 


5 


0,73 


6 


0,68 


7 


0.63 


8 


0.58 


9 


0.53 


10 


0,48 


U 


0,43 


12 


0,38 


13 


0,33 


14 


0,28 


15 


0.23 


16 


0,18 


17 


0,13 


18 


0,08 


19 


0,03 


20 


— 0.02 


21 


— 0,07 


22 


- 0,12 


23 


- 0,17 


24 


— 0,22 


25 


— 0.27 


26 


— 0.32 


27 


— 0,37 


28 


— 0,42 


29 


— 0,47 


30 


- 0.52 


31 


— 0,62 


32 


— 0,72 


33 


— 0,82 


34 


— 0,92 


35 


- 1.02 


36 


- 1,12 


37 


— 1,22 


38 


— 1.32 


39 


— 1,42 


40 


- 1.52 


41 


— 1,62 


42 


— 1,72 


43 • 


— 1,82 


44 


— 1.92 


45 


— 2,02 


46 


— 2,04 



0,98 

0,930021 

0.880160 

0,830531 

0.781241 

0.732391 

0.684077 

0.636380 

0.589394 

0,543187 

0,497832 

0.453397 

0,409939 

0,367610 

0,326163 

0,285938 

0,246874 

0,209006 

0,172361 

0.1369i33 

0,102830 

0.069982 

0.038422 

0,008164 

- 0.020794 
0,048451 

- 0,074810 

• 0,099884 

• 0,123676 

■ 0,146204 

- 0,167487 

- 0,206388 

- 0,240567 
. 0,270252 

- 0,295694 

- 0.317174 

■ 0,335008 

■ 0,335008 

- 0,361031 

- 0J69918 

- 0,376548 

- 0,581283 

- 0.384508 

- 0.386606 

- 0,387955 

- 0,388955 

- 0,389000 



1,00 

0.998796 

0,995268 

0,989538 

0,981732 

0,971970 

0,960380 

0,947081 

0,932186 

0,915821 

0,898095 

0,879128 

0,859027 

0,837897 

0,815852 

0,792991 

0,769413 

0,745221 

0,720502 

0.695357 

0,669864 

0,644120 

0,618195- 

0,592177 

0,566135 

0,540147 

0,514280 

0,488599 

0,463169 

0,438045 

0,413286 

0,365075 

0,318925 

0,275179 

0,234130 

0,196035 

0,161109 

0,161109 

0,101443 

0,076966 

0,056186 

0,039169 

0,025966 

0,016603 

0,011108 

0,009494 

0 



0 

0,047750 

0,093004 

0,135770 

0,176062 

0,2i3899 

0.249309 

0,282019 

0,ai2957 

0,341271 

0,367290 

0,391069 

0,412645 

0,432079 

0,449414 

0,464714 

0,478027 

0,489420 

0,498948 

0,506678 

0,512666 

0,516980 

0,519686 

0,520844 

0,520423 

0.518786 

0,515699 

0,511328 

0.505731 

0.498980 

0,491134 

0,472399 

0,450011 

0,424433 

OJ96102 

0,365426 

0,332789 

0,332789 

0,262988 

0,226417 

0,189076 

0,151171 

0,112873 

0,074307 

0,035575 

0 

0 



— 0,98 

— 0,930021 

— 0,889160 

— 0,830531 

— 0,781241 

— 0,732391 

— 0,684077 

— 0,636380 

— 0,589394 

— 0.543187 

— 0,497832 

— 0,453397 

— 0,409939 

— 0,367510 

— 0,326163 

— 0,285938 

— 0,246874 

— 0,209006 

— 0,172361 

— 0,1366963 

— 0,102830 

— 0,069982 

— 0,038422 

— 0,008164 
+ 0,020794 

0,048451 
0,074810 
0.099884 
0.123676 
0,146204 
0,167487 
0,206388 
0.240567 
0,270252 
0,295694 
0,317174 
0,335006 
0,335008 
0,361031 
0,369918 
0,376548 
0.381283 
0,384508 
0,386606 
0,387955 



TABLE Vlll(a).- AUXILIARY FUNCTIONS FOR INCOMPRESSIBLE JET 

FOR Baa^ 0. 



n'n 




<iF--F 


/=■'» 


a-. 

n,n 




^r-F 




0 


0,9S 


0 


1,000 


16 


— 0,6J 


- 0,0200 


0,133 


1 


0.88 


— 0,0042 


0,990 


17 


— 0,72 


— 0,0110 


0,102 


2 


0,78 


— 0.0155 


0,964 


18 


— 0,82 


0,0446 


0,07S6 


3 


0,6S 


— 0,0310 


0,922 


19 


— 0,92 


0,0803 


0,0384 


4 


0,58 


— 0,0488 


0,869 


20 


— 1,02 


0,1173 


0.0384 


5 


0,48 


— 0.0667 


0,806 


21 


— 1,12 


0,1546 


0,0259 


6 


0.38 


— 0,0835 


0,738 


22 


— 1,22 


0,1915 


0.0166 


7 


0.28 


— 0,0978 


0,666 


23 


— 1,32 


0,2272 


0,0102 


8 


0,18 


— 0,108' 


0,591 


24 


- 1,42 


0,2606 


0,0059 


9 


0,08 


— 0,114t 


0,520 


25 


- 1,52 


0,2911 


0,0031 


10 


— 0,02 


— 0,1162 


0,449 


26 


— 1,62 


0,3178 


0,0016 


11 


— 0,12 


— 0,1126 


0,382 


27 


— 1,72 


0,3398 


0,0006 


12 


— 0,22 


— 0,1037 


0,320 


28 


— 1,82 


0,3364 


0,0002 


13 


— 0,32 


— 0,0898 


0,264 


29 


- 1,92 


0,3666 


0,0001 


14 


— 0,42 


— 0,0708 


0,214 


30 


— 2,02 


0,3697 


0 


15 


— 0,52 


— 0,0474 


0,171 

: 


31 


— 2,04 


0.3890 


0 
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Figs. 1,2 




Figoxe 2.- Temperature and velocity fields in the "bo-andary layer 

of a jet at the distance x = D from the nozzle 
according to the tests of Huden. 
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Pigs. 3,4,5,6,7 
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Figure 3.- Velocity fields in the 

"boundary layer of a 
.jet at various temperatures. 
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Figure 4.- Density fields in the 

free jet at At = +60°C 
and At = -600C. 
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Figure 5.- 
of a jet. 
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Figure 5.- Field of velocity heads 

in the free .jet at At=-60°C. 
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Figure 7.- Fields of 

transverse 
velocities in the 
"boundary layer of a .jet 
at b\, = -60°C. 
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Tigs. 8,9,10,11,12,13 
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0 . .5 1.0 

Figure 8.- Velocity fields in the 

■botmiary layer of a 
high velocity Jet (Bag = l) . 
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Figure 9.- Velocity fields in the 

"boundary layer of a 
high velocity jet (BaQ = 0.5). 
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■boundary layer of a 
high velocity jet. 
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igare 10.- Density fields in the -n- tj j> t -x t. j • 

V ^ 1 -F Figure 11.- Fields of velocity heads m 
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a high velocity jet. 
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0 .5 cp 1.0 

Figure 12.- Transverse velocity 

field in a high 
velocity jet. 



0 .5 ^ 1.0 

Figure 13,- Conparison of velocity 

fields of I cooled ai^ , 
(At 60°) and II jet with high 
velocities (Bag = 1) , 
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